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APPARATUS AND METHODS FOR CHANNEL ALLOCATION 



FIELD OF THE INVENTION 

The present invention relates to apparatus and methods 
for channel allocation. 

BACKGROUND OF THE INVENTION 

State of the art methods for channel allocations are 
described in the following documents: 

"Resource allocation in wireless networks", "Scott 

Jordan, JHSN, 10 January 1995. 

Duque-Anton et al , "Static and dynamic channel 
assignment using simulated annealing", Chapter 10 in Neural 
networks in telecommunications , Ben Yuhas and Nirwan Ansari 
(Ed.), Kluwer Academic Publishers, Boston, 1994. 

"Prioritized channel borrowing without locking: A 
channel sharing strategy for cellular communications" , H . Jiang 
and S. S. Rappaport, IEEE/ ACM Transactions on Networking, Vol. 

4, No. 2, April 1996. 

The disclosures of all publications mentioned in the 
specification and of the publications cited therein are hereby 
incorporated by reference. 
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SUMMARY OF THE INVENTION 



The present invention seeks to provide improved 
apparatus and methods for channel allocation. 

There is thus provided, in accordance with a . preferred 
embodiment of the present invention, a method for utilization of 
a first plurality of channels by a second plurality of transmit- 
ters, the method including defining a third plurality of trans- 
mitter subsets such that at least one of the second plurality of 
transmitters is included in each transmitter subset, and assign- 
ing at least one channel from among the first plurality of chan- 
nels to each transmitter subset, to be shared among the transmit- 
ters in that transmitter subset, such that less than all of the 
first plurality of channels are assigned to the third plurality 
of transmitter subsets, thereby defining a reservoir of channels 
which have not been assigned to any transmitter subset, and 
sharing the channels in the reservoir of channels between all of 
the second plurality of transmitters. 

Further in accordance with a preferred embodiment of 
the present invention, a first transmitter is entitled to use a 
channel in the reservoir even if the channel is being used by a 
second transmitter, if there is no neighbor-clique which includes 
both the first and the second transmitters, wherein a neighbor- 
clique of an individual transmitter subset includes all transmit- 
ter subsets which share at least one common transmitter with the 
individual transmitter subset. 

Also provided, in accordance with another preferred 
embodiment of the present invention, is a method by which an 
individual transmitter transmits, in a situation in which a first 
plurality of channels is serving a second plurality of transmit- 
ters including the individual transmitter, the method including 
transmitting over a first channel from among the first plurality 
of channels if the transmitter belongs to a subset of transmit- 
ters which is served by the first channel and if the first chan- 
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nel is available, and otherwise, if a- reservoir of channels 
includes an available second channel, transmitting over the 
second channel, wherein the reservoir includes all channels from 
among the first plurality of channels which do not serve any 
subset of transmitters. 

Further in accordance with a preferred embodiment of 
the present invention, the channels are separated by their trans- 
mission frequency. 

Still further in accordance with a preferred embodi- 
ment of the present invention, the channels are separated by 
their transmission code. 

Additionally in accordance with a preferred embodiment 
of the present invention, the channels include CDMA (code divi- 
sion multiple access) channels. 

Further in accordance with a preferred embodiment of 
the present invention, at least some of the channels include 
wireless channels . 

Still further in accordance with a preferred embodi- 
ment of the present invention, the subset defining step also 
includes selecting, for each subset, from among the transmitters 
in the subset, a subset master to which channel allocation re- 
quests are addressed over a control channel. 

Further in accordance with a preferred embodiment of 
the present invention, the subset master is selected so as to 
maximize utilization of the control channel for communication of 
transmitters in the subset with transmitters in other subsets to 
which the subset master belongs. 

Still further in accordance with a preferred embodi- 
ment of the present invention, the transmitter in the subset 
which belongs to the largest number of other subsets is selected 

as subset master. 

Further in accordance with a preferred embodiment of 
the present invention, the method also includes releasing a 
dropout transmitter by disconnecting the dropout from the subsets 
to which the dropout belongs including informing only the masters 
of each of the subsets that the dropout has been disconnected. 
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Also provided, in accordance with another preferred 
embodiment of the present invention, is a system for utilization 
of a first plurality of channels by a second plurality of trans- 
mitters, the system including a channel assigner operative to 
assign at least one channel from among the first plurality of 
channels to each of a third plurality of transmitter subsets, 
each including at least one of the second plurality of transmit- 
ters, the channel to be shared among the transmitters in that 
transmitter subset, such that less than all of the first plurali- 
ty of channels are assigned to the third plurality of transmitter 
subsets, thereby defining a reservoir of channels which have not 
been assigned to any transmitter subset, and a channel sharer 
operative to share the channels in the reservoir of channels 
between all of the second plurality of transmitters. 

A particular advantage of the method of the present 
invention is that typically, the method is NP-complete. 



# 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention will be understood and appreciat- 
ed from the following detailed description, taken in conjunction 
with the drawings in which: 

Fig. 1 is a simplified functional block diagram of a 
channel allocation system, constructed and operative in accord- 
ance with a preferred embodiment of the present invention, which 
allocates channels to a plurality of transmitters distributed in 
space; 

Fig. 2 shows a plurality of transmitters distributed in 
Euclidian space, among obstacles, such as walls, which are also 
distributed in space; 

Fig. 3 shows the same plurality of transmitters dis- 
tributed in non-Euclidian space; 

Fig. 4 is an illustration of several subsets of 
transmitters; 

Fig. 5 shows a graph representing relationships between 
subsets I to VI of Fig. 3; 

Fig. ' 6 illustrates results of a subset graph coloring 

process ; 

Figs. 7A - 7F illustrate neighbor-cliques for each of 

subsets I - VI ; 

Fig. 8 is a table summarizing the relationships between 
the subsets and the neighbor-cliques in the example of Figs. 2 
7F; 

Fig. 9 is a table indicating which channels belong to 
each local reservoir before assignment of channels to subsets; 

Fig. 10 is a table indicating which channels belong to 
each local reservoir after channels A and B have been assigned to 
the subsets I - VI ; 

Fig. 11A is a table indicating which channels belong to 
each local reservoir after a borrowing process shown in Fig. 11B 
has occurred; 

Fig. 11B illustrates a borrowing process whereby trans- 
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mitter #1 borrows channel C from its local reservoir; 

Fig. 12A is a table indicating which channels belong to 
each local reservoir after a borrowing process shown in Fig. 12B 
has occurred; 

Fig. 12B illustrates a borrowing process whereby trans- 
mitter #21 borrows channel D from its local reservoir; 

Fig. 13A is a table indicating which channels belong to 
each local reservoir after a borrowing process shown in Fig. 13B 
has occurred; 

Fig. 13B illustrates a borrowing process whereby trans- 
mitter #15 borrows channel C from its local reservoir; 

Fig. 14A is a table indicating which channels belong to 
each local reservoir after a borrowing process shown in Fig. 14B 
has occurred; 

Fig. 14B illustrates a borrowing process whereby trans- 
mitter #17 borrows channel D from its local reservoir; 

Fig. ISA is a table indicating which channels belong to 
each local reservoir after a borrowing process shown in Fig. 15B 
has occurred; 

Fig. 15B illustrates a borrowing process whereby trans- 
mitters #1, #15, #17 and #21 have each returned the channels they 
were using and transmitter #4 has borrowed channel C from its 
local reservoir; 

Fig. 16A is a table indicating which channels belong to 
each local reservoir after a borrowing process shown in Fig, 16B 
has occurred; 

Fig. 16B illustrates a borrowing process whereby trans- 
mitter #8 has borrowed a channel from its local reservoir (trans- 
mitter #4 has not yet returned its channel; 

Fig. 17A is a table indicating which channels belong to 
each local reservoir after a borrowing process shown in Fig. 17B 
has occurred ; 

Fig. 17B illustrates a borrowing process whereby trans- 
mitters #7 and #9 each borrow a channel from their respective 
local reservoirs ; 

Fig. 18A depicts time series of transmitters #1 - #10, 
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over a relatively high channel-request-intensity time period of 
500 msec. As shown, each transmitter can either be in an "up" 
state or in a "down" state. The "up" state represents idleness 
(i.e. no channel required) whereas the "down" state represents a 
requirement for a channel; 

Fig. 18B is a "zoom" onto a 100 msec time period within 
the time period of Fig. 18A; 

Fig. 19A depicts time series of transmitters #1 - #10, 
over a relatively low channel-request-intensity time period of 
500 msec. As shown, each transmitter can either be in an "up" 
state or in a "down" state. The "up" state represents idleness 
(i.e. no channel required) whereas the "down" state represents a 
requirement for a channel ; 

Fig. 19B is a "zoom" onto a 100 msec time period within 
the time period of Fig. 19A; 

Fig. 20 is a simplified flowchart illustration of a 
preferred method of operation for the connectivity matrix 
generator of Fig. 1; 

Fig. 21 is a simplified flowchart illustration of a 
preferred method of operation for the transmitter subset 
generator of -Fig. 1; 

Fig. 22 is a simplified flowchart illustration of a 
preferred method of operation for the subset graph construction 
unit of Fig. 1; 

Fig. 23 is a simplified flowchart illustration of a 
preferred method of operation for the subset graph coloring unit 
of Fig. 1; 

Fig. 24 is a table showing input format suitable for 
either Appendix A or Appendix B; 

Fig. 25 is a simplified flowchart illustration of a 
preferred method of operation for the local reservoir management 
unit of Fig. 1; 

Fig. 26 is a simplified flowchart illustration of a 
preferred method for performing the minimal cost channel computa- 
tion step 870 in Fig. 25; 

Figs. 27A - 27B are spectral intensity diagrams of two 
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transmitters, the diagrams forming an illustrative definition of 
the terms "central channel", "adjacent channel" and "alternate 
channel" ; 

Fig. 28 is a schematic illustration of a transmitter 
which is shown communicating with other network elements in three 
different ways; 

Fig, 29 is a schematic illustration showing two trans- 
mitters in the same subset, communicating simultaneously; 

Fig- 30 is a diagram of a step-by-step process for 
subset generation; and 

Fig. 31 is an illustration of chromatic coloring 
performed over a subset graph. 

Fig. 32 is an illustration of an undirected subset 
graph with a given edge length. 

Fig. 33 is an illustration of data flow within a basic 
cell structure. 

Attached herewith are the following appendices which 
aid in the understanding and appreciation of one preferred embod- 
iment of the invention shown and described herein: 

Appendices A and B are alternative embodiments, imple- 
mented in software, of the present invention, which receive a 
transmitter connectivity matrix as input, and perform the 
following functions: transmitter subset generation, subset graph 
construction, and neighbor clique computation. 

Appendix C is a software listing of a preferred 
technique for providing a supervised optimization cycle of the 
method of Appendix A or Appendix B; 

Appendix D is an initialization file for Appendices A 

or B; 

Appendix E is an example of an output file generated by 
running Appendices A or B on an ASCII file containing the data of 
Fig. 24; and 

Appendix F is a software listing of a Matlab procedure 
performing the functions of units 60 and 70 of Fig. 1. 
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DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

A portion of the disclosure of this patent document 
contains material which is subject to copyright protection. The 
copyright owner has no objection to the facsimile reproduction by 
anyone of the patent document or the patent disclosure, as it 
appears in the Patent and Trademark Office patent file or re- 
cords, but otherwise reserves all copyright rights whatsoever. 

Fig. 1 is a simplified functional block diagram of a 
channel allocation system, constructed and operative in 
accordance with a preferred embodiment of the present invention, 
which allocates channels to a first plurality of transmitters 
distributed in space. Typically, the transmitters serve as access 
points for a larger population of workstations (not shown). 

' "The , _ term "channel" is intended to include any 
means or agency or path or line of access along which signals or 
data or information may be sent, including but not limited to 
wired channels such as TV cable and long distance underground 
cable; wireless channels such as radio channels, optical fiber 
and satellite; and any combination of the above. 

The term "connectivity" is defined as mutual 
sensitivity between any two or more transmitters. The system of 
Fig. 1 includes a connectivity matrix generator 10 operative to 
generate a matrix whose elements represent levels of 
connectivity between each pair of transmitters within the first 
plurality of transmitters. 

A transmitter subset generator 2 0 is operative to 
generate a second plurality of subsets, each subset including 
some of the transmitters. The transmitters in each subset 
typically share at least one channel which no-one outside the 
subset is entitled to use, except for transmitters outside the 
subset which can re-use the channel without substantially 
interfering with the transmitters in the subset. 

The transmitters within a subset can communicate with 
one another, e.g. when an "ordinary" transmitter asks the sub- 
set's master transmitter for a channel, by any suitable wired or 
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wireless means. Typically, one of the subset's channels may serve 
for communication within the subset. This channel is termed the 
subset's "control channel". 

It is appreciated that the transmitters in each subset 
may make any suitable use of the channel or channels shared by 
them. The channel or channels may be used solely or partly (e.g. 
by a time-sharing scheme) for some or all of the following 
functions of which (a) and (b) are control functions and (c) is a 
normal use function: 

a> Allowing a particular transmitter in the subset to 

request a channel for itself from the main reservoir, as 
described in detail below. For example, a transmitter may use 
channel B in order to request a channel for itself and in 
response may receive channel C from the main reservoir. 

It is appreciated that, alternatively, requests from a 
transmitter to the main reservoir to allocate a channel to that 
transmitter may be transmitted over a special channel which is 
used for no other purpose, i.e. is not one of the channels in the 
main reservoir. The special channel may even be over a different 
medium than the channels in the main reservoir. 

b. Allowing a workstation served by the subset to request 
access to a transmitter . For example, a workstation may use 
channel Ej'in order to announce its "desire" to send data or be 
served by the Internet. In response, one of the transmitters 
(e.g. transmitter No. 3) may use a special channel to request, 
from the main reservoir, a channel for itself. In response, the 
main reservoir may allocate channel C to transmitter No. 3 and 
the workstation may then send data to transmitter No. 3 via 
channel C. Alternatively, the request from transmitter No. 3 to 
the main reservoir may be transmitted over channel B, if a 
suitable time-sharing scheme is employed for channel B. 

c. Allowing a particular transmitter to transmit data, 
typically in order to serve a particular workstation. 

Any suitable protocol may be employed to govern 
communication between workstations and transmitters such as, for 
example, TDMA (time division multiple access) ; ring protocols 
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such as Token Ring; ALOHA ; PPMA (pre-emptive polling multiple 
access) ; GRAP (group random access polling) ; CSMA (carrier sense 
multiple access) such as CSMA-CD (CSMA-collision detection); and 
CDMA (code division multiple access) . Any suitable protocol may 
be employed to govern communication between transmitters within a 
subset , and communication between transmitters between subsets, 
such as, for example, any of the above protocols or any suitable 
combination thereof . 

The process implemented in Appendix A is preferably 
performed before the operational stage is active. 

Referring back to Fig. 1, a subset graph construction 
unit 3 0 is operative to construct a graph to represent 
connectivity of the subsets, where subsets are considered 
connected if their respective closest member transmitters are 
relatively close to one another. 

A neighbor clique computation unit 34 is operative to 
compute the set of all neighbor cliques in the subset graph 
constructed by unit 30. A preferred method of operation for unit 
34, including an example of suitable source code, is described in 
pages 59 - 69 of Lau, H. T. , Algorithms on Graphs , TAB Profes- 
sional and Reference Books, Blue Ridge Summit, PA, USA, 1989. 

A graph coloring unit 40 is operative to assign 
channels to the subsets by coloring the subset graph generated by 
unit 30, the color of each subset (node) representing a channel 
assigned thereto. Channel assignment by graph coloring allows a 
minimal number of channels to be used, while ensuring that use of 
a channel by each subset does not substantially interfere with 
the transmitters in any other subset. A functional block diagram 
of a preferred embodiment of graph coloring unit 4 0 is 
illustrated in Fig. 23. A preferred embodiment of a graph 
coloring unit, including an example of suitable source code, is 
described in Chapter 4 of Lau, H. T., Algorithms on Graphs . TAB 
Professional and Reference Books, Blue Ridge Summit, PA, USA, 
1989. Algorithms on Graphs also provides code for performing 
subset graph coloring. 

A reservoir manager 50 is operative to manage a main 
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reservoir including all channels which have not been assigned to 
specific subsets. The reservoir manager 50 is a functional unit 
which, in practice, is typically performed by individual 
transmitters. The individual transmitters performing the 
reservoir management function may be predetermined but preferably, 
they are elected. In the illustrated embodiment, a subset master 
assignment subunit 60 is operative to assign or elect transmit- 
ters to perform the reservoir management function, typically one 
master transmitter per subset. Preferably, all of the master 
transmitters are of equal status and there is no "master-master 
transmitter". In certain applications, however, it may be desira- 
ble to assign or predetermine a "master-master transmitter", 
which may or may not comprise one of the master transmitters. 

Typically, the master transmitter of each subset is 
selected as the transmitter in the subset which belongs to the 
largest number of other subsets, i.e. the transmitter having the 
largest neighbor clique. It is appreciated that if the master 
transmitter of a subset A is also a member of a subset B then 
preferably, the transmitters within subset A can utilize subset 
A's control channel in order to transmit messages to members of 
subset B, via subset A's control channel, the master transmitter, 
and the control channel of subset B. This is typically the case 
even if the master of subset A is not the master of subset B. 
Preferably, this utilization is only permitted if the control 
channel is idle, i.e. requesting channels from the master subset 
takes precedence over general communication with other subsets 

via the master subset. 

A particular advantage of selecting as subset master, 
the transmitter with the largest neighbor clique is that general- 
ly, this selection maximizes the connectivity between subsets 
by' maximizing utilization of the control channel with consequent 
efficient use of the non-control channels allocated to that 
subset. 

Each manager transmitter, serving a particular subset, 
includes a local reservoir management unit 70, which is a 
functional unit operative to manage the main reservoir from the 
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point of view of that subset- The main reservoir, as managed by a 
particular subset manager on behalf of a particular subset, is 
also termed herein the "local reservoir" of that subset, as 
described in detail below. It is appreciated that the subset 
manager itself typically sometimes requires a channel and in this 
case, typically, the subset manager requests a channel from 
itself similarly to the way other transmitters in the subset 
request a channel from the manager. 

The term "assign" is used to refer to the process of 
providing a subset with a channel which "belongs" to that subset 
and is not merely "borrowed" by the subset. The term "allocate" 
is used to refer to the process of "loaning" a channel to an 
individual transmitter from the main reservoir. 

Fig. 2 shows a plurality of transmitters 110 
distributed in Euclidian space, among obstacles 120, such as 
walls, which are also distributed in space. 

Fig. 3 shows the same plurality of transmitters 
distributed in non-Euclidian space. It is appreciated that the 
metric used to define distances between transmitters is typically 
a non-Euclidian metric. For example, if a non-Euclidian space 
such as a Labochevsky or Minkowski metric is employed, the dis- 
tance between transmitters is proportional to the signal strength 
propagation (typically measured in electromagnetic dB) existing 
between those transmitters. The Minkowski metric is described in 
Space Time Motion , Moyseyev, R. , Chapter 2 , Dover Books, 1942. 
Uses of the Minkowski metric are described in P. G. Bergmann, 
Introduction to the Theory of Relativity , Dover, Prentice-Hall, 
1942. A suitable Minkowski metric for the illustrated embodiment 
of the present invention is described in the following website of 
the Department of Physics at Virginia State University: 
http : //astsun . astro . Virginia , edu/ " eww6n/math/MinkowskiSpace . html . 

The plurality of transmitters is typically such that 
there is substantially no cross-talk i.e. the distance between 
each transmitter and itself is infinite which, under the 
Labochevsky metric is equivalent to a zero distance. 

It is appreciated that the transmitters of Fig. 3 may 



13 



be regarded as a complete graph wherein each transmitter is a 
node and the weight of the edge connecting each two nodes is the 
distance between them. The distance may, for example, be a Labo- 
chevsky metric distance, expressed in dB. If two transmitters are 
positioned such that they cannot receive from one another, the 
distance between them is typically considered infinity. 

Transmitter subsets are now defined such that at least 
one of the 21 transmitters (in the present example) is included 
in each transmitter subset and, preferably, each transmitter is 
included in at least one subset. 

Fig. 4 shows 6 subsets, each including 5 or 6 
transmitters although it is appreciated that more generally, any 
number of transmitters , including only one transmitter, may be 
included in a particular subset. 

To generate the subsets of Fig. 4, the following 
criteria may be employed: 

a< There is a predetermined maximum number of transmitters 

per subset and a predetermined minimum number of transmitters per 
subset. For example, the rule may be that each subset must 
include between 5 and 6 transmitters, as in Fig. 4. 

b. For- each transmitter, there is a maximum number of 

subsets to which that transmitter is entitled to belong. In the 
illustrated example, each transmitter is entitled to belong to 
only up to 2 subsets. However, more generally, some transmitters 
may be entitled to belong to only, say, 3 subsets, whereas 
another transmitter may be entitled to belong to, for example, 
only one subset. 

In order for a transmitter to belong to n subsets, the 
transmitter typically has the ability to operate in n different 
channels simultaneously, e.g. typically includes n separate and 
independent transmitter units e.g. n different radios. For 
example, in order for a radio transmitter to belong to 2 subsets, 
the radio transmitter typically includes two separate and 
independent radio units. This is one reason for establishing a 
maximum number of subsets to which each transmitter is entitled 
to belong. However, in some applications, it is desirable for the 
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maximum number of subsets for a particular transmitter, to be 
less than the number of channels through which that transmitter 
is capable of transmitting. 

c. The number of transmitters which belong to only a 

single subset should be minimized. 

<3- The number of subsets which are connected should be 

minimized. Two subsets are "connected" if the distance between 
them is below a predetermined connectivity threshold such as, in 
the illustrated example, -91 dB. 

The "distance" between two subsets with no common 
members is defined as the minimum distance between any 
transmitter in the first subset and any transmitter in the second 
subset. The distance between subsets having a common member is 
zero and therefore, subsets having a common member are always 
considered connected. 

e. The number of neighbor-cliques to which each subset 

belongs, should be minimized. Each subset S has a "neighbor- 
clique" which includes all subsets which are connected to S. 
Therefore, the number of neighbor-cliques to which each subset S 
belongs is equal to the number of subsets which are connected to 
S. 

Any suitable technique may be employed to generate 
subsets which answer to the above criteria (a) to (e) , such as 
techniques based on combinatorial algebra or neural networks or 
artificial life methods or Monte Carlo methods. Publications 
describing artificial life methods include: 

D. Whitley, "An executable model of a simple genetic 
algorithm". in Foundations of genetic algorithms — 2, L. D. 
Whitley (Ed*)/ Morgan Kaufmann, San Mateo, CA, USA. 

M. D. Vose, "Modelling simple genetic algorithms", in 
Foundations of genetic algorithms — 2, L. D. Whitley (Ed.), 
Morgan Kaufmann, San Mateo, CA, USA. 

Next , a graph is generated to represent relationships 
between the subsets. Typically, each node in the graph represents 
a subset and a pair of nodes is adjacent (i.e. there is an edge 
between them) if the corresponding subsets are connected. Fig. 5 
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shows a graph representing relationships between subsets I to VI 
of Fig. 3. 

The available channels are color-coded and the nodes 
(subsets) in the graph are now colored, using a minimum number of 
colors. In other words, the color of a node identifies a channel 
belonging to the subset corresponding to that node. It is 
appreciated that if more than one channels has been assigned to a 
particular subset, that subset's node will have more than one 
color (i.e. the length of the color vector of that node will be 
more than one) . Any suitable method may be used to color the 
channels such as chromatic coloring, described in Chapter 4 of H. 
T. Lau, Algorithms on graphs , TAB Professional and Reference 
Books, Blue Ridge Summit, PA, USA. 

The channels which are not assigned to any subset are 
regarded as belonging to a main reservoir which is used commonly 
by all the transmitters in all the subsets. As described in more 
detail below, channels in the main reservoir can sometimes be 
used simultaneously by more than one transmitter. This situation 
is termed herein "channel reuse 11 . 

Channel re-use is typically facilitated by deriving 
local reservoirs from the main reservoir. One local reservoir is 
typically derived for each transmitter subset. The channels 
within a local reservoir at a given time t are always a subset of 
the channels included in the main reservoir at time t. A channel 
is included simultaneously in two local reservoirs I and II iff 
that channel can be used simultaneously by any transmitter in the 
transmitter subset I corresponding to reservoir I and by any of 
the transmitters in the transmitter subset II. corresponding to 
reservoir II. More generally, a channel is included simultane- 
ously in n local reservoirs iff that channel can be used simulta- 
neously by one member each of the n transmitter subsets of the n 
reservoirs respectively without substantial interference. 

Typically, two transmitters are considered to be 
operating "without substantial interference" if the value of a 
predetermined cost function of the transmission quality between 
the two transmitters is below a threshold value. For example, the 



16 




cost function may comprise the BER (bit error rate) of the 
channel between two transmitters . 

Initially, each of the local reservoirs includes all of 
the transmitters in the main reservoir. However, it is 
appreciated that subsequently, the local reservoirs each shrink 
and expand back again differently, depending on each subset 1 s 
relationship with other subsets and depending on the extent to 
which each subset needs channels. 

An example of the development of 6 local reservoirs 
from a main reservoir serving 6 subsets of transmitters is now 
described, based on the example of Figs . 2 - 6 . 

Figs. 7A - 7F illustrate neighbor-cliques 200, 210, 
220, 230, 240 and 250 for each of the subsets I - VI 
respectively. These cliques are defined by neighbor clique compu- 
tation unit 34 of Fig. 1. Each clique is identified by a dotted 
line surrounding all of its members. 

Fig. 8 is a table summarizing the relationships between 
the subsets and the neighbor-cliques in the example of Figs. 2 
7F . "x" indicates connectivity between two subsets . 

Fig. 9 is a table indicating which channels belong to 
each local reservoir before assignment of channels to subsets as 
shown in Fig. 6. As illustrated, each of the local reservoirs 
includes all of the channels A - G and therefore, of course, the 
main reservoir also includes all of the channels A - G. 

Fig, 10 is a table indicating which channels belong to 
each local reservoir after channels A and B have been assigned to 
the subsets I - VI as shown in Fig. 6. As illustrated, each of 
the local reservoirs now includes only channels C - G and there- 
fore, of course, the main reservoir also includes only channels C 
- G. It is appreciated that at this stage, all six local reser- 
voirs are identical . 

Fig- 11A is a table indicating which channels belong to 
each local reservoir after transmitter #1 has borrowed channel C 
from its local reservoir, as shown in Fig. 11B. Since transmitter 
#1 (cross-hatched) belongs to subset II, as shown in Fig. 11B, 
and since subset II is connected to only one other subset (subset 
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I) , channel C is removed only from the local reservoirs of sub- 
sets I and II and remains present in all other local reservoirs, 
i.e. the local reservoirs of subsets III - VI. 

Fig. 12A is a table indicating which channels belong to 
each local reservoir after (as shown in Fig. 12B) transmitter #21 
has borrowed channel D from its local reservoir, and before 
channel C has been returned by transmitter #1. Since transmitter 
#21 also belongs to subset II, as shown in Fig. 12B and since 
subset II is connected to only one other subset (subset I) , 
channel D, like channel C, is also removed only from the local 
reservoirs of subsets I and II and remains present in all other 
local reservoirs, i.e. the local reservoirs of subsets III - VI. 

Fig. 13A is a table . indicating which channels belong to 
each local reservoir after (as shown in Fig. 13B) transmitter #15 
has borrowed channel C from its local reservoir, even though 
channel C has not been returned by transmitter #1. This is 
possible because transmitters #1 and #15 belong to subsets II 
and VI respectively, which are not connected. Since transmitter 
#15 belongs to subset VI, as shown in Fig. 13B and since subset 
VI is connected to only one other subset (subset IV) , channel C 
is removed from the local reservoirs of subsets VI and IV and 
remains present in the local reservoirs of subsets III and V. 

Fig. 14 A is a table indicating which channels belong to 
each local reservoir after (as shown in Fig. 14B) transmitter #17 
has borrowed channel D from its local reservoir. This is 
possible, even though channel D is being used by transmitter #21 
because the subsets VI and II to which transmitters #17 and #21 
belong are not connected. Since transmitter #17 belongs to subset 
VI, as shown in Fig. 14B and since subset VI is connected to only 
one other subset (subset IV) , channel D is removed only from the 
local reservoirs of subsets IV and VI and remains present in the 
local reservoirs of subsets III and V. 

Fig. 15A is a table indicating which channels belong to 
each local reservoir after (as shown in Fig. 15B) transmitters 
#1, #15, #17 and #21 have each returned the channels they were 
using and transmitter #4 has borrowed channel C from its local 
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reservoir. Transmitter #4 belongs to subsets I and III, as shown 
in Fig. 15B. Since subset I is connected to subsets II, III and 
IV and subset III is connected to subsets I and V, channel C is 
removed from the local reservoirs of subsets I - V and remains 
present only in the local reservoir of subset VI. 

Fig. 16A is a table indicating which channels belong to 
each local reservoir after (as shown in Fig, 16B) , transmitter #8 
has borrowed a channel from its local reservoir (transmitter #4 
has not yet returned its channel) . 

Fig. 17A is a table indicating which channels belong to 
each local reservoir after (as shown in Fig. 17B) , transmitters 
#7 and #9 each borrowed a channel from their respective local 
reservoirs (transmitters #4 and #8 have not yet returned their 
channels) . 

Fig. 18A depicts time series of transmitters #1 - #10, 
over a relatively high channel-request-intensity time period of 
500 msec. As shown, each transmitter can either be in an "up" 
state or in a "down" state. The "up" state represents idleness 
(i.e. no channel required) whereas the "down" state represents a 
requirement for a channel. 

Fig. 18B is a "zoom" onto a 100 msec time period within 
the time period of Fig. 18A. 

Fig. 19A depicts time series of transmitters #1 - #10, 
over a relatively low channel-request-intensity time period of 
500 msec. As shown, each transmitter can either be in an "up" 
state or in a "down" state. The "up" state represents idleness 
(i.e. no channel required) whereas the "down" state represents a 
requirement for a channel. 

Fig. 19B is a "zoom" onto a 100 msec time period within 
the time period of Fig. 19A. 

Fig. 20 is a simplified flowchart illustration of a 
preferred method of operation for the connectivity matrix genera- 
tor 10 of Fig. 1. The input to the method of Fig. 20 is an 
indication of the locations of the transmitters in space, e.g. 
using a conventional Euclidian metric system. Each transmitter 
performs carrier measurement for all other transmitters for a 
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specified time period such as approximately 2 seconds for a 
typical 8-level FSK (frequency shift keying) modulation system 
for the ISM band under FCC (Federal Communications Committee) 
Part 15 regulations. 

Preferably, in addition to the carrier measurements, 
the correlation detection rate (CDR) of the transmitters is 
measured, e.g. by giving a FOM (figure of merit) for existence of 
clock within the carrier. Unit 330 couples between conjugate 
information, i.e. information regarding the same transmitter, 
arriving from units 310 and 3 20. The output of unit 3 0 is fed to 
a connectivity matrix computation procedure 334 which comprises 
looped blocks 340, 350 and 360. 

The connectivity matrix generator 10 is operative to 
generate a matrix whose elements represent levels of connectivi- 
ty between each pair of transmitters in the system. In the illus- 
trated example in which there are 21 transmitters, the connectiv- 
ity matrix is a 21 x 21 symmetrical matrix whose diagonal is 
zeros . 

Fig. 21 is a simplified self-explanatory flowchart 
illustration of a preferred method of operation for the transmit- 
ter subset generator 20 of Fig. 1. In step 410, criterion (a) is 
the criterion (a) described above with reference to Fig. 4. 

In step 430, "nearest" typically is determined by 
computing, for each unassigned transmitter, the shortest distance 
between it and any assigned transmitter from among the various 
assigned transmitters, and selecting as "nearest" the unassigned 
transmitter having the smallest "shortest distance". The metric 
used in this step is not Euclidian but rather a metric of the 
connectivity matrix e.g. Minkowski. 

In step 440, criteria (a) - (e) are the criteria (a) 
(e) described above with reference to Fig. 4. 

Optionally, step 460 is performed. In this step, the 
assignment of transmitters accomplished in steps 400 - 450 is 
optimized using conventional methods such as evolutionary genetic 
algorithmic methods such as those described in the Vose publica- 
tion referred to herein. 
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The output of the method of Fig. 21 is typically a 
transmitter subset vector such as that appearing in Appendix E. 

Fig. 22 is a simplified flowchart illustration of a 
preferred method of operation for the subset graph construction 
unit 30 of Fig. 1. The method of Fig. 22 receives a transmitter 
subset vector such as that generated by the method of Fig. 21, 
and treats each subset as a node in a graph. The method adds 
edges to connect each pair of subsets which have a common 
transmitter, thereby to generate a transmitter subset graph such 
as the graph of Fig. 5. 

Fig. 23 is a simplified flowchart illustration of a 
preferred method of operation for the subset graph coloring unit 
4 0 of Fig. 1. The input to the method of Fig. 2 3 is the subset 
graph generated by the method of Fig. 22. In step 610, an 
"available" color is a color not yet assigned to any node. The 
output of Fig. 23 is a colored subset graph. An example of a 
colored subset graph is illustrated in Fig. 6. 

Fig- 25 is a simplified state machine illustration of 
the local reservoir management unit 70 of Fig. 1. Each 
transmitter has two states: idle, and "needing a channel". 

In step 830, the term "subset master" refers to the 
master of the subset to which the transmitter belongs or, if the 
transmitter belongs to more than one subset, to any of the 
masters of any of the subsets to which the transmitter belongs. 

Fig. 2 6 is a preferred method for performing the 
minimal cost channel computation step 870 in Fig. 25. 

Figs. 27A - 27B form an illustrative definition of the 
terms "central channel", "adjacent channel" and "alternate chan- 
nel". Figs. 27A and 27B are spectral intensity diagrams of first 
and second transmitters respectively. As shown, the first trans- 
mitter's spectrum (Fig. 27A) includes a central lobe 1010, pri- 
mary side lobes 1020 and 1030 and secondary side lobes 1040 and 
1050. The second transmitter's spectrum (Fig. 27B) includes a 
central lobe 1060, primary side lobes 1070 and 1080 and secondary 
side lobes 109 0 and 1100. The central lobe of each transmitter 
occupies a channel which is termed the "central channel" of that 



21 




transmitter. When the specification speaks of a transmitter 
"borrowing a channel", what is meant is that the borrowed channel 
becomes that transmitter's central channel. 

The two primary side lobes of each transmitter occupy 
two respective channels which are termed the two "adjacent chan- 
nels" of the transmitter. The two secondary side lobes of each 
transmitter occupy two respective channels which are termed the 
two "alternate channels" of the transmitter. Therefore, if the 
central channel of one transmitter A is the adjacent channel or 
even the alternate channel of another transmitter B, then there 
is interference between transmitters A and B as shown in Figs. 
27A - B. 

A "zero cost" borrowing, from the point of view of 
transmitter A relative to transmitter B, refers to a situation 
wherein B borrows, a channel but neither of that channel's adja- 
cent channels or alternate channels overlap with any of A's 
central, adjacent or alternate channels. 

A non-zero cost borrowing, from the point of view of 
transmitter A relative to transmitter B, refers to a situation 
wherein B borrows a channel but at least one of that channel's 
adjacent channels or alternate channels overlap with at least one 
of A's central, adjacent or alternate channels. 

Typically, transmitter A is only permitted to borrow a 
channel if, for all other transmitters T it is true that the 
interference introduced by the borrowing process is no greater 
than the ratio between the amplitudes of T's central channel and 
either of T's adjacent channels. 

The level of interference caused to transmitter A by 
transmitter B is termed "adjacent channel interference level" if 
the ratio between A's central channel's amplitude and B's central 
channel's amplitude is equal to the ratio between the amplitudes 
of A's central and adjacent channels. 

Similarly, the level of interference caused to trans- 
mitter A by transmitter B is termed "alternate channel interfer- 
ence level" if the ratio between A's central channel's amplitude 
and B's central channel's amplitude is equal to the ratio between 
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the amplitudes of A's central and alternate channels. 

Typically, if two channels are in the same reservoir, 
there is no interference between them. 

Fig. 28 is a schematic illustration of a transmitter 
1500 which is shown communicating with other network elements in 
three different ways, as indicated by outgoing communication 
channels A, B and C (dotted lines), and incoming communication 
channels D, E and F (solid lines) . 

As shown, transmitter 1500 transmits information, via 
channel A, to an individual network element 1510' from among a 
plurality of network elements 1510 which are within a wireless 
communication envelope around transmitter 1500, defining those 
network elements with which transmitter 1500 is capable of commu- 
nicating. Transmitter 1500 transmits information, via channel B, 
to another network element 1520, which is illustrated as being in 
the same subset 153 0 as transmitter 1500 although this need not 
necessarily be the case. Channel B is a channel allocated to 
transmitter 1500 in accordance with a preferred embodiment of the 
present invention. Transmitter 1500 additionally transmits infor- 
mation, via channel C and via a subset master 1540 common to both 
subsets 1530 and 1550, to a transmitter 1560 in subset 1550. 

Transmitter 1500 receives information, via channel D, 
from an individual network element 1510'' from among the plurali- 
ty of network elements 1510 which are within the wireless commu- 
nication envelope around transmitter 1500. Transmitter 1500 
receives information, via channel E, from another network element 
1520, which is illustrated as being in the same subset 1530 as 
transmitter 1500 although this need not necessarily be the case. 
Channel B is a channel allocated to transmitter 1500 in accord- 
ance with a preferred embodiment of the present invention. Trans- 
mitter 1500 additionally receives information, via channel F and 
via a subset master 154 0 common to both subsets 153 0 and 1550, 
from a transmitter in subset 1550 such as transmitter 1560 or any 
other transmitter in subset 1550. 

Fig. 29 is a schematic illustration showing two trans- 
mitters 1500 and 1540 in the same subset 1530, communicating 
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simultaneously (with transmitters 1540 and 1520 respectively) . 
This is possible because one of the transmitters 1500 is using 
the channel allocated to subset 153 0 whereas the oth-^-r 
transmitter (subset master 1540) is using the control channel of 
subset 1530 . 
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A preferred process for generating subsets is now 
described . 

Given a plurality of Transmitters (Access Points) 
randomly distributed over a bounded domain, some of which 
may be connected to wired networks, the Transmitters may 
possibly interfere with each other, and cause 
interference in the network. 

The following description defines the sapient 
network segmentation and access engines involved in 
presenting a segmented topology that provides maximum 
network efficiency and capacity, while preserving maximal 
frequency reuse and robustness i.e. builds a network 
presenting a strong connectivity between the Transmitters 
with maximal network data flow. 

The following approach clusters the Transmitters 
into logical Subsets, whose physical characteristic is 
that of being a complete graph (where the Transmitters 
are the graph vertices).'. Following this Subset 
segmentation of the transmitter plurality, a Subset graph 
emerges. In order to enhance performance, the emerged 
graph should fulfill some requirements (discussed below). 

Assuming a limited and insufficient number of 
frequencies available, graph coloring is used for 
assigning a minimal needed frequency regime. to the Subset 
graph. This process is called FCA (Fixed Channel 
Assignment). In operational terms, a local frequency bank 
is defined (the frequency reservoir set), and Subsets 
shall borrow frequencies from this set upon channel 
demand . 

In a multi-Subset structure, with a number of 
inherent topological characteristics, mobile stations 
present free dynamics in regard to Subset coverage and 
each mobile station preferably associates with at least 
one Subset from the cluster. Access is performed via a 
slotted-aloha mechanism of channel request acquisition 
followed by R-TDMA data sequencing. Every transmitted 



J:\Richard\IIca\Insert 7/12/98 



25 



data entity accepted by a specific transmitter in the 
Subset, shall be analyzed for its destination resolving. 
Subsequent decisions shall forward this data entity over 
the Subset (which is at this point is assumed to be 
wired) to its destination transmitter (the transmitter 
responsible for the final access step to the 
destination ) . 

Initial and Boundary Conditions include : 
a transmitter set, {AP} ; 

an inter transmitter interference map, when dB 
is the mapping matrix, and dB H is the absolute value of 

the interference ratio for transmitter ( I ) and 
transmitter ( j ) ; 

a defined threshold for the relative 
interference ration given in absolute value dB ; 

a minimal number of Transmitters on a Subset, 

a maximal Number of Transmitters on a Subset, 

u max ' 

a' maximal number of radio interfaces on every 
transmitter: p; 

a finite set of. discrete frequencies, {/} ; and 
a number of fixed channels for a Subset, k. 

Topology Generation 

The following described model, maps the transmitter 
set into a Subset set {R} , while preserving: 

1) boundary conditions: d, e and f 

2) minimization of Transmitters that belong to a single 
Subset, this can be expressed by minimization of 

1l*ll. aj , -a^ i * / - \ : 1 

3) Maximization of Subsets that are not interfering, 
thus frequency re-use can be implemented. The non- 
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interference is defined as 
<P = {AP}\aP eR f rsR J J = \ i ..,R,j= 1,..,/?,/* j 

4) Minimization of cliques that a Subset belongs to. 

When a neighbor clique is defined as 

{C,. } = {AP}\APeR f r>R j J = \ 7 ..,R,i* j 

The rank of a given Subset topology is defined as 
min(Cos/ la ...,Cosr / ) when Cost, is defined as Cost; - Cost c +■ co • Cost^ , 



Cosi c = min 



Vi=i J 



Cost > ■= xx\\n $ , and d is the number of 



different constellations (possible combinatory), and co is 
a weighting factor. 

Subset generation is assumed to follow each transmitter 
birth (potential energy transmission), while transmitter 
death will cause no Subset changes until some pre-defined 
cost function violation (this cost function can be 
throughput criterion, connectivity loss, extreme external 
interference, etc.). 

A channel assignment engine running over the system 
performs assignment (graph coloring) over the Subset 
topology, thus., as graph cliques contain fewer Subsets, 
the number of fixed frequencies (vertex colors) required 
to color the graph, becomes lower. 

Fig. 30, and the following description, is a step by step 
model for Subset generation under the initial conditions, 
and fulfilling the above listed requirements. 

Initially an AP k is chosen that is located at "mass 
center" of the transmitter cluster, i.e., satisfies the 
following condition, 

^-^^minjX/^f, When N ~^ AP ^' 
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The first Subset is constructed by aggregating 
Transmitters to the first assigned using the following 
manner, 

7?, = {AP t } + {APX:;VB, <^,V/ = s + \..N. dB b SdB™ 

The remaining Subsets are build in two iterated 
steps, step 1 choose an transmitter that is the nearest 
transmitter to all ready existing Subsets, step 2 
aggregates Transmitters that would be connecting- 
Transmitters between previous Subsets and the current 
one. And step 3, choose additional free Transmitters 
until the defined bound. 

Step 1: Choose an AP k that is located at "mass center" of 
the transmitter cluster, i.e., satisfies the following 
condition, 

]T. dB kj = min|^ y t/5 iy j When j runs over the Transmitters in 

existing Subsets, and AP k is chosen from the free. 
Transmitters. 

Step 2: Choose {AP}™ in such that each APj belongs to an 
existing Subset, and APj preserves condition f, and there 
is no other transmitter in any of the existing Subsets 
that its dB to APk is smaller than the dB of Apj to Apk. 

The Subset is completed to a max f rom the free 
Transmitters set via steps 1 and 2. 

HCA 

Assuming no spreading, the system is required to 
solve a channel assignment problem. In general channel 
assignment rating is weighted by its spectral span (that 
is preferably the minimum needed) and the blocking 
probability (blocking is defined as the condition when a 
channel is needed but not available). 
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The channel assignment process can be split into two 
different operational regimes (thus this is a Hybrid 
process ) : 

Fixed Channel Assignment ( FCA ) . that assigns a 
given number of channels to each Subset k, and by so 
doing ensures Subset traffic connectivity. 

Dynamic Channel Assignment (DCA) - that serves 
as a channel reservoir to Subset channel requests. As 
Subset channel request rate is spatio-temporal (a 
function of time and domain, i.e., Subset) different, 
reservoirs will evolve during an operational regime. 

E£A 

Chromatic-Coloring is performed over the Subset 
graph by assigning k colors to each Subset, and by 
mapping those colors into the channel set, see Fig. 31 
for illustration of coloring. After so doing, a new 
channel set is defined (the set of non-assigned channels) 
and name it as the channel reservoir set. 

A coloring of an undirected graph G is an assignment of 
colors to nodes of G such that no adjacent nodes of G 
have the same color. The chromatic number of G is the 
minimum number of colors needed to cover G. 

IWjhllWolh^'^fe) When C «( G n) is the chromatic number 
for G R . 

The coloring is performed by a simple implicit 
enumeration tree search method. initially, node 1 is 
assigned color 1, and the remaining nodes are colored 
sequentially so that node I is colored with the lowest- 
numbered color which has not been used to color any nodes 
adjacent to node 1. Let p be the number of colors 
required by this feasible coloring. Attempt to generate 
a feasible coloring using q<p colors. To accomplish this 
all nodes colored with p must be recolored. Thus, a 
backtrack step can be taken up to node u, where node u+I 
is the lowest index assigned color p. Attempt to color 
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node u with its smallest feasible alternative color 
greater than its current color. If there is no such 
alternative color which is smaller than p, then backtrack 
to node u-2. Otherwise, recolor node u and proceed 
forward, sequentially recoloring all nodes u+2,., 
with the smallest feasible color until either node n is 
colored or some node v is reached which requires color p. 
In the former case, an improved coloring using q colors 
has been found; in this case backtrack and attempt to 
find a better coloring using less than g colors. In the 
latter case, backtrack from node v and proceed forward as 
before. The algorithm terminates when backtrack reaches 
node 1 . 



|{/U 



DCA 

As defined above, after the FCA procedure there are 
available channels in the reservoir. {/} res is 

assigned to each Subset in G R . Each Subset borrows 
channels from this set upon internal request. The 
criteria for letting a channel from the reservoir is 
defined as follows: 

A graph of nearest-neighbors for a given Subset request 
- G nn is defined thus, 

G nn = {K,£}when V = {/?,,...,/?„} and £ = such that 

J^ntRj *0|/* j 

The available channels from{/}^to a given Subset is 
defined as {/}„ = {/}„ - f/}^ when {/}^ = {/}„„. -{/}„„, . 

The blocking probability for the proposed channel 
borrowing based DCA is now defined, and thus 
demonstrating the Subset topology generation model 
stipulations. 

{/}* = {/L-y{/L^/U when vM«*-. is the unified set 

of all channels used in the cliques or Subset R and 
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{f} coil ±s the set of collided channels (discussed later). 
In other words ^{f} cij(Jue - {/} se!f "^{/Ij 1 represents neighbor 



requests, and 



ue i ~ ^j a ^ - From these notations, it is 



clear that {/} diqUL , depends upon traffic in the neighbor 
cliques. It is appreciated the number of channels 

preferably needed to construct jj{/} re Jj in such a way that 
no blocking shall appear in the system. 

After FCA the R Subset needs as peak demand a R - 1 
channels, thus in order to eliminate the blocking 



9~ i| 



probability 



«.-i = K/L|-2>.-K/} 



\ KJ 'coll* 



and 



that 



From this, minimizing a R or 
\c\ will lead to the minimization of ||{/} r< J| t . , remembering 
that ||{/} r J| ..is a vector on R . 



Assuming j 



{/} 



is a given constant, then system 
level optimum is given by minimization of ^ {/}, 

; = l 

Thus the blocking probability P B ( ) R can be given by 

/ % . If/Li . 1 

P B [) R = mm 



jj v-' / res R j 



Note : the stiffness of the 

r res * II J 

probability decay is a function of the traffic demand and 
the access mechanism. 

Access 

As described for Subset-generation, the Transmitters 

are addressed as vertices of a wired graph, and. are 

clustered into highly-connected Subsets. Logically each 
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Subset serves as a network segment, thus networking-wise, 
each Subset becomes a wireless cell. Stations (mobile 
clients) wandering under the coverage of a specific 
Subset receive and transmit information (in a packet 
based fashion) via the Subset "gateways" - the 
Transmitters. It is of no significance to the station 
which transmitter collects its information, and how this 
information passes on the network. 

A two layered protocol is now described: 

the in-Subset protocol, also referred to as the 
radio MAC protocol, that is responsible for coordination 
over the radio transactions and for message delivery; and 

the Subset-protocol that is responsible for 
Subset -transmitter control , frequency assignment and 
inter-Subset connectivity . 

Following is a general description, and the joint 
mathematical model, of the protocols. 

The described access method fulfills two major 
tasks: over the coordination, and over the air high- 
capacity data transaction. Each of the tasks is bounded 
to some boundary conditions: the coordination function of 
the MAC is bounded to perform transmission request 
acquisitions and Transmitters-on-Subset timing 

synchronization. The transaction part of the MAC is 
referred to as LLC (Link Layer Control) and its purpose 
is to handle channel requests (from frequency assignment 
point-of -view ) , and to handle data traffic to and from 
the radio channel. 

Every transmitter on the Subset separately performs 
a Slotted-Aloha based acquisition of access requests and 
then handles transmissions under an R-TDM regime. In 
order to perform the R-TDM section of the protocol, the 
DCA algorithms are activated by the requesting 
transmitter (these requests are sent over the Subset 
media, to the Subset master . which is one of the 
Transmitters ) . 

Periodically, the transmitter accepts a permission 
to generate a slotted-aloha competition window. Using a 
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dedicated channel, it sends a CCLR broadcast packet 
(competition-start announcement) over the radio. In 
response, all active stations not engaged with some 
transmitter compete in this process by sending a random 
access based RTS (Request To Send) control frame. The 
transmitter stores all access requests in a request- 
buffer (RTS queue). Using the specified queue as a FIFO, 
it polls the awaiting stations for data transmission. In 
case stations are addressed from the on-Subset side, 
reverse polling is used. 

As slotted aloha mechanisms are very sensitive to 
the correctness of the randomization seed, the 
transmitter should calculate the estimated number of 
active stations in the Subset at the next aloha 
competition window by advising with rest of the 
Transmitters on the same Subset, and then include this 
knowledge to the potential competitors over via the CCLR 
control frame. 

The second layer protocol performs three main tasks: 
the token control passing for the slotted-aloha 
competition timing for the transmitter Subset population, 

serves as the back-bone for channel assignment 
coordination: requests and releases of frequency basic- 
bands; and 

handles and transfers data traffic between the 
Transmitters in the Subset neighbor clique. 

This is a token based protocol, and two types of 
tokens may exist over the same Subset. The Subset 
coordination token whose purpose is to coordinate 
slotted-aloha competition timing. On receiving such a 
token, the transmitter sends back on the Subset a channel 
request from the Subset master (a different or the same 
transmitter) in case an bver-the-radio transmission is 
needed (RTS queue is not empty, or there is available 
data for its radio stations). The master preferably 
responds immediately by either granting or denying the 
channel (at this point all neighbor-Subsets population is 
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announced). It is noted that this token rolls just 
within its home Subset, and is transparent to the 
neighbor Subsets in the neighbor clique. The second token 
type rolls in the neighbor clique, i.e., circles over all 
neighboring Subsets. It serves like a token-Subset access 
grant, that upon receiving it , a non-idle transmitter is 
allowed to transmit its data packets (up to some pre- 
defined burst length given in time units). 

Destination address resolution is performed in the 
following manner. Upon receiving a data packet (invariant 
to its source location), the transmitter compares it to 
the Subset current population, i.e., every transmitter is 
assumed to have complete information regarding current 
address-queue status of all Transmitters on the same 
Subset. In case the destination address is located on 
the Subset (no information on neighbor Subsets), the 
transmitter forwards this packet to the home-transmitter, 
otherwise this packet is addressed as non-resolved, and 
is being transferred to the neighboring Subset in such a 
way this packet will circle the complete neighbor clique. 

Each transmitter performs IP routing from the 
Subset-topology viewpoint. Every data packet transmitted 
over the Subset topology, is encapsulated with its 
initial radio transmitter (home transmitter) IP address. 
The destination transmitter Subset ( transmitter Subset 
that the destination station is located behind, could be 
a wired LAN or a mobile station) strips this 
encapsulation, and maps the source transmitter IP to the 
MAC packet address. 

In order to avoid packet loopings and Broadcast 
storms, a minimal spanning tree self configuration is 
implemented. 

As system conditions change as a function of time- 
dependent interference and varying traffic regimes, the 
spanning tree graph is to be re-computed periodically. 
The objective is defined as following: Consider an 
undirected graph G with given edge length (in our case 
traffic amount and its typical BER). The objective is to 
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find a spanning tree in G such that the sum of the edge 
lengths in the tree is minimum, see Fig. 32. 

The following process achieves this objective over G 
(Subset graph). Initially, set T is empty. Edges are 
considered for inclusion in T in the non-decreasing order 
of their lengths (costs). An edge is included in T if it 
does not form a cycle with the edges already in r. A 
minimum spanning tree is formed when n-1 edges are 
included in. T . in the implementation, the edges are 
partially sorted with their smallest edge at the root of 
a heap structure (a binary tree in which the weight of 
every node is not greater than the weights of its sons). 
The running time of the above process is 0(m log m) , 
where m is the number of edges in the graph. 

The basic cell structure is described as follows and 
is illustrated in Fig. 33. The arrows indicate data 
flow. It is assumed that the queue behavior of an 
investigated transmitter is governed by six stochastic 
processes : 

A. Queue addition via in-Subset air protocol 
(regardless of destination location;. 

B. Queue addition via the Subset protocol, from 
another transmitter located on the same Subset, i.e., the 
destination address is in the investigated transmitter 
queue; 

C. Queue addition via the Subset protocol, from 
another transmitter not located on the same Subset, i.e., 
the destination address is in the investigated 
transmitter queue; 

D. Queue subtraction via in-Subset protocol, i.e., 
destination located behind the investigated transmitter, 
and is reflected in the transmitter's queue; 

E. Queue subtraction via Subset protocol, i.e., 
destination located behind an transmitter located on the 
same Subset ; and 

F. Queue subtraction via Subset protocol, i.e., 
destination located behind an transmitter located on the 
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same Subset, but its destination transmitter is not a 
Subset member. 

A hidden Markov model for queue dynamics is assumed. 
The system level model is constructed from an assembly of 
those queues. Queue dynamics is determined via stochastic 
equations, computed in event-step-jump time . The 
derivation of these equations and their mutual coupling 
is shown as follows. 

{Q} 0 =*{P) = {dQ}=>{QL 
[dQ] =>{(&}„ 



When {Q} u = {Q}-*+{<*Q} 



1 ) Construction of A process 

The basic entry probability for a single slotted-aloha 

r 



1 



1 



1- — When n K 



slot by a single station is given by 

is the estimated number of competing stations, and n a is 
the actual competing station number. Given a probability 
distribution V(n) specifying the probability that any two 
stations in the Subset are at. some dB threshold 
difference, i.e., partial interference over the same 
slot, the probability of slot drawing success for a u 
station collision over the same slot could be evaluated 

■r*(n) When k = — -; Thus, the complete slot 

v J 2 



— f 1- — 



access probability under BER will be given by: 
p = S-Y — {\-fer)\ 1- — [l-fer) -r*(n); When S is the number 
of slots in the slotted-aloha window . 
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Number of active stations is given by: n a = y^Q RTS . ; 

When y^2/?r^ ^ s "the number of stations in the RTS queues of 
the Transmitters in the Subset, 

The estimated number of active stations is preferably 
presented via some algorithmic form. As the estimation 
error increases, the access latency through the aloha 
mechanism increases non-linearly . 

Given the. probability of a station entering the 
transmission request queue it is preferably translated 

into an actual number of packets: A0 R7 - S = min(l, ; As 
time steps for the probability estimation is A/ CC£J? which 
the aloha token passing duration, and it is very short, 
AQ rts cbu not be immediately translated into the packet- 
queue. Thus, virtual arrival times are calculated and 
assigned. Let us define an arrival time vector 

Q--{^l&Z when '°='"~^; p ' is the 

probability of granting a frequency via the DCA process 
and t prev ±s the arrival time of the previous packet 
(initial time is the floating time t). 

2) Construction of the B and C processes. Define P DP to be 
the probability of a packet originating in a given 
Subset, to be destined to a station in the same Subset. 
The probability that the destination address is located 
in the investigated transmitter's queue is given by 



the assumption in this expression is that the 



destination addresses of packets, are distributed 

(r ) 

uniformly over the address domain when tr " is the total 
number of stations in the Subset (active and idle). 
Thus, the number of packets the investigated transmitter 
is preferably able to handle is given by 
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^=(>t)l (p - &); when 



is the investigated 



transmitter. 

The C process . queue addition due to packet incoming 
from remote Subsets in the same neighbor clique. The 
probability for the investigated transmitter to handle 



such a packet is given by 



O 3 



Qc 



5C3-1 



(0 



3C3-1 

When y\R {i) 
(=1 



is the 



number of Subsets in the neighbor clique, not including 
the investigated transmitter's Subset. The assumption 
made here is that all Subsets operate under the same 
operational regime. 

Combining the two processes (B an C), provides the 
number of packets generated within the complete neighbor 
clique that are to be handled by the investigated 
transmitter in the At CCLR interval- Given a Subset-token to 
a transmitter, it can transmit over the Subset for a 

Br 



defined time period given by A/ 7 



Rate,, 



When Br is the 



burst length in minimal-packets. Thus the number of 
packets to be handled by the investigated transmitter is 



given by &Q B + C = 



P Q. +P Qc 

ica 

Z* (t) 



At, 



CCLR 



A/ T 



Rate, 



When 



A/, 



At 



■ ■ Rates 



77? 



IS 



the number of packets transmitted during the A/ CC£J? duration 
(this could be a packet fraction). 

3) Construction of the E and F processes. The probability 
that the investigated transmitter needs to transmit to 
the neighbor clique (not to its own air) is given by 

PE,F-^PoR^-^-PoR\Qirs and after simplifying this 



expression 



)ing 



this 



probability expression into number of packet is given by 
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o J 



1 QrTs[ ^DR '[ n (R f ) 



+ 1-1 



Thus, number of packets 



transmitted by the investigated transmitter to the 
neighbor clique (not over its own Subset) is given by 



0 J 



1-O j \ P ■ 

L ^RTS [ 1 DR 



1 



n 



+ 1-1 



N'HQrts 

k-\ 



Mr 



•Rate D 



; When IV is the 



number of Transmitters over the entire neighbor clique 



oca 



and is given by ^ = 2^ - 



k = \ 
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The above description included herein assumes ring 
protocols for communication between transmitters within a subset 
and recursive slotted ALOHA for communication between 
transmitters between subsets. 

A particular advantage of a preferred embodiment of the 
present invention is that the suitability of the apparatus is not 
limited to communication systems having a particular transmitter 
distribution in interference space. Typically, conventional 
systems are suitable for a particular, symmetric transmitter 
topology. According to a preferred embodiment of the present 
invention, the apparatus is suitable for substantially any 
transmitter topology and in particular for applications in which 
transmitters are mobile and therefore, transmitter topology 
varies with time. 

Another particular advantage of a preferred embodiment 
of the present invention is that transmitters belonging to more 
than one subset are particularly suited to serve as relays, i.e. 
to relay information between the subsets to which they belong. 

The maximum number of transmitters per subset is typi- 
cally determined by the subset transmission rate. For example, if 
the speed of the subset's medium is 100 mbps and each transmitter 
is capable of transmitting at most 10 mbps, then no more than 10 
transmitters are typically allowed to belong to each subset in 
order to prevent waste of channels. 

A preferred software implementation of the present 
invention is set forth in Appendices A - C which contain computer 
listings of a software embodiment of the present invention. 

Appendices A and B are alternative embodiments, 
implemented in software, of the present invention, which receive 
a transmitter connectivity matrix, such as that generated by unit 
10 of Fig. 1, as input, and perform the functions of units 20, 
30, 34 and 40 of Fig. 1. 

Appendix C is a software listing of a preferred 
technique for providing a supervised optimization cycle of the 
method of Appendix A or Appendix B. 

Preferred methods for providing supervised optimization 
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cycles are described in the above-referenced Vose publication. 

Appendix D is an initialization file for Appendices A 

or B; 

Appendix E is an example of an output file generated by 
running Appendices A or B on an ASCII file containing the data of 
Fig. 24; 

Appendix F is a software listing of a Matlab procedure 
performing the functions of units 60 and 70 of Fig. 1. The 
procedure runs on Matlab (Matrix Laboratory) for Windows, Version 
5.0 and up, commercially available from The MathWorks Inc., 
Cochituate Place, 24 Prime Park Way, Natick, MA 01760, USA. The 
procedure receives as input an output of Appendices A or B, such 
as the example output of Appendix E. 

Fig. 24 is a table showing input format suitable for 
either Appendix A or Appendix B. The contents of the table of 
Fig. 24 is an example of input which suitably represents the 
example of Figs. 2 - 19B in the present specification. 

A preferred method for utilizing the computer listings 
of Appendices A - C is as follows: 

a. Using a PC 486 equipped with Delphi Pascal, generate a 
unit and key in the contents of Appendix A and either B or C, 
under the appropriate headings (Unit, interface, uses, const, 
etc.) as shown in Appendices A - C. 

b. Compile and run the software using the "main" program 
of Delphi Pascal. 

It is appreciated that the software components of the 
present invention may, if desired, be implemented in ROM (read- 
only memory) form. The software components may, generally, be 
implemented in hardware, if desired, using conventional tech- 
niques. 

It is appreciated that the particular embodiment de- 
scribed in the Appendices is intended only to provide an extreme- 
ly detailed disclosure of the present invention and is not in- 
tended to be limiting. 

It is appreciated that various features of the inven- 
tion which are, for clarity, described in the contexts of sepa- 
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rate embodiments may also be provided in combination in a single 
embodiment. Conversely, various features of the invention which 
are, for brevity, described in the context of a single embodiment 
may also be provided separately or in any suitable subcombina- 
tion. 

It will be appreciated by persons skilled in the art 
that the present invention is not limited to what has been par- 
ticularly shown and described hereinabove. Rather, the scope of 
the present invention is defined only by the claims that follow: 
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CLAIMS 



1. A method for utilization of a first plurality of 
channels by a second plurality of transmitters, the method 
comprising: 

defining a third plurality of transmitter subsets such 
that at least one of the second plurality of transmitters is 
included in each transmitter subset; 

assigning at least one channel from among the first 
plurality of channels to each transmitter subset, to be shared 
among the transmitters in that transmitter subset, such that less 
than all of the first plurality of channels are assigned to the 
third plurality of transmitter subsets, thereby defining a 
reservoir of channels which have not been assigned to any 
transmitter subset; and 

sharing the channels in the reservoir of channels 
between all of the second plurality of transmitters. 

2. A method according to claim 1 wherein a first 
transmitter is entitled to use a channel in the reservoir even if 
the channel is being used by a second transmitter, if there is no 
neighbor-clique which includes both the first and the second 
transmitters, wherein a neighbor-clique of an individual 
transmitter subset comprises all transmitter subsets which share 
at least one common transmitter with the individual transmitter 
subset . 

3. A method by which an individual transmitter transmits, 
in a situation in which a first plurality of channels is serving 
a second plurality of transmitters including the individual 
transmitter, the method comprising: 

if the transmitter belongs to a subset of transmitters 
which is served by a first channel from among the first plurality 
of channels, and if the first channel is available, transmitting 
over the first channel; 
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otherwise, if a reservoir of channels includes an 
available second channel, transmitting over the second channel, 
wherein the reservoir includes all channels from among the first 
plurality of channels which do not serve any subset of 
transmitters . 

4. a method according to claim 1 wherein said channels are 
separated by their transmission frequency. 

5. a method according to claim 1 wherein said channels are 
separated by their transmission code. 

6. A method according to claim 5 wherein said channels 
comprise CDMA (code division multiple access) channels. 

7. a method according to claim 1 wherein at least some of 
said channels comprise wireless channels. 

8. A method according to claim 1 wherein the subset 
defining step also comprises selecting, for each subset, from 
among the transmitters in the subset, a subset master to which 
channel allocation requests are addressed over a control channel. 

9. a method according to claim 8 wherein the subset master 
is selected so as to maximize utilization of said control channel 
for communication of transmitters in the subset with transmitters 
in other subsets to which the subset master belongs. 

10. A method according to claim 8 wherein the transmitter 
in the subset which belongs to the largest number of other sub- 
sets is selected as subset master. 

11. A method according to claim 1 and also comprising 
releasing a dropout transmitter by disconnecting the dropout from 
the subsets to which the dropout belongs including informing only 
the masters of each of the subsets that the dropout has been 
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disconnected . 



12. A system for utilization of a first plurality of chan- 
nels by a second plurality of transmitters, the system compris- 
ing : 

a channel assigner operative to assign at least one 
channel from among the first plurality of channels to each of a 
third plurality of transmitter subsets, each including at least 
one of the second plurality of transmitters, the channel to be 
shared among the transmitters in that transmitter subset, such 
that less than all of the first plurality of channels are as- 
signed to the third plurality of transmitter subsets, thereby 
defining a reservoir of channels which have not been assigned to 
any transmitter subset; and 

a channel sharer operative to share the channels in the 
reservoir of channels between all of the second plurality of 
transmitters. 

13. A system according to claim 12 and substantially as 
shown and described above. 



14. A system according to claim 12 and substantially as 
illustrated in any of the drawings. 

15. A method according to any of the preceding claims 1 
11 and substantially as shown and described above. 

16. A method according to any of the preceding claims 1 - 
11 and substantially as illustrated in any of the drawings. 



gor the Applicant, 



Sanford T. Colb & Co, 
C:29465 
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APPENDIX A 



unit Rinq_ce 




on routine; 




interface 



uses 

Windows, Messages, Sysucils , Classes, Graphics, Controls, Forms, 
Dial c g s , 

StcCtrls, ExtCtrls; 



TFcrml = class (TForm) 
5u:tcnl : TButtcn; 
3ulICh2 : TButrcn; 
Panell : TPanei; 
Edit!: TEdit;- 
Buttcn3: TBut:on; 

procedure But toniClick ( Sender : TObject) ; 

procedure Sutton2Click ( Sender : TObject) ; 

procedure But tonSCIick ( Sender : TObject) ; 
private 

{ Private declarations } 
public 

■[ Public declarations } 



const 



MaxNumof A? = 500; 
^axNumOf Rings - 50; 



type 



end; 



r crml : T: or ml ; 



{I/O related variables} 



X_File 
Mass File 



string; 
s t ring; 
string ; 
string; 
string; 
string ; 
string; 



R_Fiie 
V_File 
A_File 
inif lie 



WriteOut 
Ini_f ile 
Out File 



Text File; 
TextFile; 



{ Common variables} 



Numo fAP 
MaxAP_onRing 
Ma :< Ri n g_ f o rAP 
Sensitivity 
Rings 



Integer ; 
Word; 



Word; 
Word; 
Word; 



Ring 



array [ 1 MaxNumOf Rings , 1 MaxNumo fAP ] of Integer 



procedure Initiation; 
procedure Main_Loop; 
procedure Genera te_Rings ; 
procedure Output; 



implementation 



{ SR ' . OFM } 




procedure rr-.a i n_lccp 
begin 

initiate; 
Ger:6raie_Rif:CS ; 
Output ; 

end; 

procedure Initiation; 

var 

i, j : Word; 
begin 

Randomize ; 

inifile := 'd: \sim\MAC-Sirn\SGHz\rinc . ini * ; 

AssicnFiie (Ini_Fiie, inifile) ; 
Reset { ±ni_Fiie) ; 

Reacln ( Ini_File, Numof AP) ; 

Readin ( Ir.i_File, Sensitivity) ; 

Readin { Ini_FiIe, MaxA?_onRinc ) ; 

Readin (Ini_Fiie / MaxRing_f orAP ) ; 

Readin {Ini_Fiie) ; 

Readin {Ini_Fiie, K_Fiie) ; 

Readin ( Ini_File,' WriteOut ) ; 
CioseFiie ( Ini File); 



{read the K Mac} 

Ass icnFile (Ini_File, K_File) ; 

Reset ( Ini_Fiie) ; 

for i := I to Numof AP do 
bee in 

for j := i to Numof AP do read ( Ini_File / K_Map[i,jj); 
readin <Ini_File) ; 
'end; 

CioseFiie (Ini File); 



end ; 



procedure Generate_rings ; 
var 

Temp, i :Word; 

Distance : array [ 1 . . MaxAP, 1 . . MaxAP] of Word; 
Temp_Distance, Pointr :array [1.. MaxAP] of Word; 
Numb f Rings : Word; 



begin 

NumOf Rings : =0; 

{ CALCULATE THE DISTANCE MATRIX FROM AP LOCATION MATRIX - THE EUCLIDEAN- 
CASE} 

for AP_a:=l to Numof AP do 
for AP b:=l to NumOfAP do 

Distance [ AP a, AP b ]: =Sqrt ( Sqr ( Field [AP_a , 1 ]- Field [APJo, 1 ] ) 



_ -r Sqr ( Field [.J^^l 2] -FieId[A?_b, 2} ) ) ; 

for i:=I to NumofA? dc 

for AP_a:=I to NumcfAP-I dc 

fcr A?Jc: =A?_a-^l to NumOfA? dp 

if DisLance[i, A?_a] > Distance [i, A?_b] then 
bee in 

Temp :=Di stance [i, A?_b] ; 
Distancefi, AF_b j : -Distance [ i , A?_a ] ; 
Di stance i i , AP_a]:=Temp; 
end; 

for A?_a:=i tc NumcfAP do 

for A?_b:=I tc NumOfA? dp 



•FIND THE A? THAT IS THE MASS CENTER OF THE CLOUD } 

for A?_a : =1 to NumOfA? do Temp_Dist ance [ Ar_a ] : =C ; 
for A?_a:=i to NumOfA? do 
for A?_b:=i tc NumOfA? do 

Temp_Qi stance [A?__a ] : =Temp_Dis tance [A?_a j -f Di s t ance '[ A?_a , AP_b] 

?emp:=l; 

for' AP _ a : -2 to NumCfAP do 

if Temp_Distance [ A?_a] < Temp_Di s t ance [ Temp ] then Temp : =A?_a ; 

{ GENERATE FIRDT RING} 
Pointr [ i] : -1; 
Ring [ 1 , 1] : =Temp; 
Inc ( RingsOnAP [Temp] ) ; 
.for AF__a:=l to MaxAP_onRing- 1 do 
begin 

fncfPomcr [I] ) ; 
Ring [I, 



procedure Output ; 
var 

i, j : Word; 
begin 

AssignFile (Out_File, WriteOut) ; 
Rewrite (Out_File) ; 

fcr i : = 1 to Rings do 
begin 

write (Out_File, 'Ring' , i: 3, * : * ) ; 
for j := 1 to NumofAP do 

write (Out_Fiie / K_Map [ Numof AP+ 1 -i , j ] ) ; 
writeln (Out__File ) ; 
end; 

CioseFiie (Out_File) ; 

end; 



procedure TForml . But toniClick { Sender : TObject); 
begin 

Pane 11. visible := true; 
Main_Loop; 



? = r.eil . Visible a ise; 

er.c; 



crccecure TFcriuI . But tcn2Ciick ( Sender : TObiect ) 
Hale; 

end; 

procedure TTorrril . Bun ton2Ciick ( Sender : 70b j ec: ) 
becin 

inifiie := edit i . text ; 
Buzz o p. I . V isicle := true; 
5ucicn2 . visible := true; 
Icitl .Visible := false; 
3uzzc-2 .Visible : - false; 

end ; 



end. 



APPENDIX 




uses 

Windows, Messages, SysUtils, Classes, Graphics, Cor.- rc Is , Forms, 
Dialogs , 

StdCiris, ExtCtris; 



cons: 



MsxNu.TiCfA? = 200; 

KaxNumofRings = BO; 

MaxNuTtcfRincs3 = 210; 

K MaxNumof A? = 50; 



type 

TFormi = class (7 Form) 

Bution: : TButton; 

Buziop.2 : T3ut ton ; 

Paneil: TPanei; 

Editl: TEdit; 

Buicon5 : TButton; 

procedure ButtcnlCiick {Sender : 

procedure Button2Ciick( Sender : 

procedure But tor.3C lick (Sender : 
private 

{ Private declarations } 
public 

{ Public declarations } 
end ; 



TCbiect ) ; 
TObj ect ) ; 
TObject ) ; 



var 

For ml 



TFormi; 



{I/C related variables} 



K_Fiie 
Mass_Fiie 
R_Fiie 
V_File 
A_Fiie 
inif ile 
WriteOut 
Ini_f ile 
Out File 



string; 
string; 
string; 
strings- 
string; 
string; 
string; 
TextFile, 
TextFile, 



{ Common variables} 



Numof AP 
MaxA?_onRing 
MaxRing_f orAP 
Sensitivity 
Rings 

MinA?_onRing 
R 

Termination 

R_?oint 

First_In_Ring 

Free_Pointr 

Treshold 

Option 



Word; 
Word; 
Word; 
Word; 

array { 1 .. MaxNumof Rings , 
Words- 
Word; 
Boolean; 
array [ 1 . 
array [ 1 , 
Words- 
Single ; 
Word; 



. H MaxNumof A?] of Word 



. MaxNumof Rings ] of Word; 
.MaxNumof Rings ] of Word; 



Di s t ance_?c 

?.incs_Cr.A? 

X _:-;£? 

C:=:ance, d 
Graph 



# 



Were; 



array [ I . . MaxMumo f A? 1 of Hyte; 

array [I . . MaxNurnof A?, 1 . . Ma xNumcf API of Wore; 

array [ I . . MaxMumo f A?, I . . MaxNuiacf A? j of Single; 

array [ I . . MaxNumo f A? ] cf Were; 

array [ 1 .. MaxNumo f Rings , 1 MaxNumo f Rings j cf Word; 



procedure Gener a t e_lcca t ions ; 

procedure Initiation; 

procedure Gener at e_rir.es ; 

p rocecure Generat e_?.i nc_G raph ; 

procedure Output; 

procedure Main_Lcop; 

implementation 



{ 3R - . DFM} 



procedure ma in_looe ; 
begin 

initiation; 
Generat e_r ings ; 
Generate_Rinc_Graph ; 
Output; 

end ; 



procedure Generat e_Locat ions ; 
var 

i, j : Word; 

begin 

Randomize ; 

AssignFile (Ini_File, K_File) ; 
Rewrite ( Ini_Fiie) ; 

Case Option of 

1: for i := 1 to NumofAP do 

writeln (Ini_Fiie, Random ( 300 ) +1 , * ' , Random ( 300 )+ 1 ) 

2: for to Trunc ( Sqrt (NumofAP) + 1 ) do 

for j:=i to Trunc (Sqrt (NumofAP) 4-1) do 

writein (Ini_Fiie, 10+i*10, f ',I0+j*10); 
end; {of Case} 

CloseFile (Ini_File) ; 

end; 



procedure Initiation; 



Generate a working matrix f o m 



procedure WcrkAre^^i generate a working matrix for^^Fe distances} 
var 

i, j : Wore; 
begin 

fcr i:=I to NumofA? cc 

for to NumcfA? do d { i , j j : -Distance [ i , j ] ; 

end ; 



procedure FreeAF; {Find the free-A?s} 
var 

i, j, k:wcrc; 
K_Nct_Free : Boolean; 

begin 

for k:=I to NumOfAP do 
begin 

K_Not_Free : =False; 
for i : =1 to R do 

for i:=i to R_?oint[i] do 

If Rings [i,jj=:< then K_Nct_Free : =True ; 



if net (?;_Not_Free) then 
begin 

Inc (Free_?oini:r) ; 

rree_Lis: [free_?om:r] :=k; 
end; 



{check if all APs are assigned} 
'if Free_?ointr = 0 then Termination : =True 
else Terminat ion : =Faise ; 

end; 

var {of ring-generation procedure} 
i , j , Min_Indx f Min_Incx_to : Word; 
k,i,m :Wcrd; 
Left: Word; 
Min : Single; - 

Temp, Temp__indx : array [ 1 . . MaxNumof AP ] of Word; 



begin 

Min_Indx : =65000 ; 
Wor kArea ; 



{find Cluster center} 

for i:=I to NumofAP do Temp[i]:=0; 

for i : =1 to NumofAP do 

for j:=l to NumofAP do Temp [ i ] : =Temp [ i ] +Trunc ( 1+d [ i , j ] ) ; 

Min:=65000; 

for i:=l to NumofAP do 
if Temp[i]<Min then 
begin 

Min: =Temp [i] ; 
Min_Indx : =i ; 

end; 



Inc(R); 

Inc (R_?oint [R] ) ; 

Rings [R, R_Pcint [R] ] : =Min_Indx; 



First__In_^^[Min_Indx] :=1; 

{ update number of rings per A?} 
Inc ( Rings_CnA?[Min_Indxj ") ; 

■ {add the res- of the A?s tc the ring} 
for j:=i tc MsxA?_cnRinc-l dc 
begin 

Min:=65000; 

fcr i:=I to NumofA? do 

if (c [Rings [ R , 1 ] , i ] <Mi n } and 
{d[Rings [ R , 1 ] , i ] > 0 } then 
begin 

'Min:-d[Rings [ R, 1] , i j ; 
Min_Indx:=i; 

end; 

Inc (R_?oint [Rl } ; 

Rings ( R, R__?cint [R] ] :=Min_Indx; 
d [ Rings [R, 1 ] , Hin_Incx] :=G; 

{update number of rings per' A?} 
Inc {Rings_OnA? [Min_Indx] ) ; 
end; 

{GENERATE ALL OTHER RINGS } 

{find the nearest Free -A? to one of the rings APs } 
Termination : =Ealse ; ■ 

While not {Termination) do' 
begin 

Wor kArea ; 

Eree_Pointr : =0 ; 
Free A?; 

{For all Free APs find the A? that has the closest NOT-Free 

A?} 

Min:=65000; 

fcr i:=I to Free_Pcintr do 
for j:=l to NumofAP do 

if (d [Free_List [i] , j ] <Min) and (d [ Free_List [i] , j ] >0) and 
(First_In_Ring [i] =0) and { Rings_OnA? [ j ] >C ) then 
begin 

Min : =d [ Free_List [ i ] , j ] ; 
Min_Indx: =i; 
end; 

if Min<65000 then 
begin 

Inc (R) ; 

Inc (R_Point [R] ) ; 

Rings [ R, R_Point [R] ] :=Free_List [Min_Indx] ; 
d [Rings [R, 1] , Free_List [Min_Indx] ] : -0; . 
First_In_Ring [Free_List [Min_Intix] ] :=1; 

{ update number of rings per AP } 
Inc (Rings_OnAP [ Free_List [Min_Indx] ] ) ; 
end; 




connected 



su:s "he 
ring 

the 



and 



iperforn^e first stace cf ceneraticn; cWcc=:e al- 

A?s up to min treshoic} 
for j : = 1 to MinA?_cnRing do 
begin 

Min:=65G00; 

for i:=i to NumcfA? do 

if (d[Rings[R, 1] , i]<tfin) and ( Rings_OnA? [ i 3 >0 ) 

(Rincs_CnA? [i ] <MaxRing_f orA? ; and 
(dfRincs !* R, I j , i J >0) then 
bee in 

Min:=dfRings [ R, I ] , i ] ; 
Min_Incx:=i ; 

end ; 

if Min<65000 then 
begin 

Inc (R_Point [R] ) ; 

Rings [ R, R_?cint [ R] ] :=Min_mc:<; 

d[ Rings [R, 1 j , Min_Indx] :=0; 

{update number of rings per A?} 
Inc (Rings_CnA? [Min_Indx j ) ; 
end; 

end; 

(perform the second stage cf Generation: co-Iecate ail free 
A?s up to max treshcld } 

{calculate treshoic: the maximal distance to a connected A?} 
{this parameter should be changed tc a better criterion to ■ 

non-Euclidian spase. the best will be if the ap in a given 

will be at the most eaual distance from one to the other. 
There is a cossibiiitv that a cetrain ap will not be able to 
connect, thus it should generate a seperate ring, anc skipp 

first stage of the connected ap acquisition} 
Treshoic: =0; 

for j:=2 to R Point [R] do t . . , 

if distance [Rinas [R, 1], Rings [R, j]]>Treshoia unen 
Treshold:=cistance [Rings [R, i] , Rings [R, j ] 1 ; 

While ( j<=MaxAP_onRing-MinA?_onRing-l) or (Termination) do- 
begin 

Min:=65000; 

for i:=l to NumofA? do 

if ( d [ Rings [ R, i j , i] <Min) and ( Rings_OnAP [ l ] -0 ) 

( d [ Rings [ R, i] , i] >0) and . 
( d [ Rings [ R, 1] , i] <Treshold) then 
bea in 

Min : =d f Rings [ R, 1] , i] ; 
Min_Indx:=i; 

end; 

if Min<65000 then 
becin 

Inc (R Point [Rj ) ; 



Rings [ R, R_?cint [Rl ] :=MJ^^hdx; 

d [Rings [R, Ij .Minjncx] 

{update number of rings per A?} 
Inc {Rings_CnA? [Min_Indxj ) ; 
end; 
Inc(j); 

end; 

if Not (Termination; then 
begin 

Free_?ointr : =0; 

FreeA?; 
end; 

end; (cf While} 



procedure Generate_Ring_Graph; 
var 

i, j, k, Count :Word; 

Graph_Temp : array [ 1 . . MaxNumof Rings 3 ] of Word; 
begi n 

for i:=i to R do for j:=i to R do Graph [ i, j ]: =0 ; 
for i:=I to MaxNumof Rings 3 do Graph_Temp f i } : =0 ; 
Ccunc:=G; 

for k:=l to NumOfAP do 

{if Rings_OnA?[k] >1 then} 
begin 

for i : =1 to R do 

for j:=l to R_Point[i] do 
if" Rings [i, j ]=k then 
begin 

Inc (Count ) ; 
Graph__Temp [Count ] 
end; 

if Counc>=2 then 

for j:=2 to Count do 
begin 

Graph [Graph_Temp [ 1] , Graph__Temp [ j ] ] r =1 ; 
Graph [Graph_Temp [j ] , Graph_Temp [ 1 ] ] :=I; 
end; 
Count : =0 ; 
end; 

end; 

procedure Output; 
var 

i, j : Word; 
begin 

AssignFile (Out_File, WriteOut) ; 
Rewrite (Out_File ) ; 

Writein (Out_Fiie , 'APs to Ring Mapping, and Ring connectivity 
list ' ) ; 

Writein (Out File ) ; 



for i : = 1 to R do 
begin 



write {^Friie, ' Ring ' , i : 3, 1 : 1 ) ; 

for j := I to R_?ointii] do write <Out_Fiie, R ir.cs [ i , j ] : 4 ) ; 

{write out the rinc connectivity information) 

for i:=i to (Ma>:A?^cnRing-R_?cint ( i ; ) do wr i t e (Cut_r i le , ' 

write (Out_Fiie, * ' ; ; 

for to R co if Graph : i , j ]= i then Write {Cut_r iie , j : 4 ) ; 

writein (Out_Fiie) ; 
end; 

: 

Writein (Cut_Fiie) ; Writein (Out_Fiie) ; 

Writein {Out_Fiie, 'Ring Graph Connectivity Matrix'}; 
Writein {Out_Fiie) ; } 
for i : =1 to R do 
begin 

write {Out_Fiie, ' Ring * , i : 2, ' : ' ) ; 

for j:=I to R do if Graph[i,jj=i then Wr it e <Out_r i ie , 1:4); 
writein (Gut_Fiie) ; 
end; } 

Close File (Cut_File) ; 

end; 



procedure TForn\i . Butt oniCiick ( Sender : TObject ) ; 
begin 

Panel i . visible := true; 
Mci"_Locp; 

Panel! . Visible := false; 

end; 

procedure TForml . Sue ton2Ciick ( Sender : TObject) ; 
begin 

Halt ; 

end; 

procedure TForml . But tcn3Ci ick ( Sender : TOb j ect } ; 
begin 

inifile := edit 1 . text ; 
Button! . Visible := true; 
3uttcn2 . Visible := true; 
Edit! .Visible : = false; 
5uttcn3 . Visible : = false; 

end; 



end . 



APPENDIX C 



prccra m G.A; 
uses crt, 



ccns: 



PocMaxSize 
GeneMaxSize 
Mut a t ionSeeci 
Thr es ho ic2 Bread 
N umO fGeneracicns 
AiieiDimension 

>iut at ionSchemaBin 

Wr i - e Frequency 
Gene rat ion 
PcpSize 



ailel 



GeneSize 
PopCount r 
scanning . 

AlielCoup.tr 
Crossover Point r 



the 



Next PopCount r 



Tot Fitness 
PopFi tness 



Of fsprificLiSL 
information . 

McdelGene 
MutationSchema 

crossing 

Gene rat ionOut File 
f 



Maximal number of genes in the system 
Maximal number of Aiieis in every gene 
Mutation rate for gene level 
Selection Strength 

Maximal number of generations in a search 
Numerical diversity of each 
Aliei [0, Ailel Dimension j 

Time bin for mutation schema calculations 
in terms of number of generations. 
The frequency of output schedule. 
Generation count r. 

A matrix containing the genes and their 

information . 
Number of aiieis. 

General purpose ccuntr for population 

Counter for the Gene level scanning. 

The location of the gene breaking point for 

purpose of crossover. 

The number of surviving individuals after 

selection stage. 

Mean population fitness. 

An array containing the fitnesses of each- 
in the population. 

An array containing the offsprings 
A golden gene . 

Low- path filter over mutation to fitness 

on an allel basis scale. 
Global log file. 
Output file. 



type 



PopMaxS i ze 

GeneMaxSi ze 

MutationSeed 

AllelMutat ion Seed 
Thres ho id2 Breed 
Allel Dimension 
• NumOf Generations 
MutationSchema Bin 

Wri te Frequency 



= 4 00; 
30; 
10; 

8; 
5; 
20; 
- 2000; 
1; 

= 100; 



field = array [1 .. PopMaxSize] of single; 

GenePop = array [ 1 PopMaxSize, 1 .. GeneMaxSize] of integer, 
TargetGene = array [ 1 .. GeneMaxSize] of word; 



var 



Gene re^^ph were! 

PcpSize : word ; 

Gene Size : word ; 

PopCcuntr : word 

A 1 1 e I C o u r. t r : word 

Crossover Point r : word 

Next PopCount r : word 

Tot Fitness : single 

PcpFitness : Field ; 

PcpList : GenePop 

Of f sorincList : Gene Pop 

ModeiGene . : TarcetGene; 

MutaticnScherrta : TarcetGene; 

GenerationOutFile :. text ; 

f : text ; 



Procedure NuilArrays; 
begin 

for PcpCountr : = 1 to PopSize do 
begin 

PopFitnessCPopCountr] := 0; 

for AiIeiCcup.tr := 1 tc GeneSize do 

Cf f springList [ PopCcur.ir , AlieiCcuntr] : = 0; 

end; 

For AlieiCcuntr : = 1 to GeneSize do 

MutaticnScherr:a [AlieiCcuntr] : = 0 ; 

end; 

{ } 

Procedure Genera teModelGene ; 
begin 

for AllelCountr := 1 to GeneSize do 

ModeiGene [ A.iieiCount r ] Random (Al lelDimens ion + 1); 

end; 

{ } 

Procedure Initiate 10; 
begin 

Ass ign ( Genera t ionOut Fi le , ' c : \algrthms\ga\ga . out ' ) ; 

Rewrite ( Generat ionOut File ) ; 

ClrScr; 

end; 

{ _ } 

Procedure WriteOut; 
begin 

writein ( Generat ionOut File ) ; 

writeln ( Generat ionOut File , 'Generation: Generation, 

Mean Fitness: ' , Tot Fi tness : 7 : 5 ) ; 
writeln (GenerationOutFile) ; 
writeln (GenerationOutFile) ; 

for PopCountr : = 1 to Popsize do 
begin 

for AllelCountr := 1 to GeneSize do 

write (GenerationOutFile, Pop-List [ Popcountr, 



^ llelCCUntr?:5 Uitein(Ger,erarionOutFiIe, ■ 
• ,?ccFitness [PopCountr) :7:5) ; 



end; 

writeln (GeneraticnOu^ile] 



jrlteir: (General ionOut t He) ; 

end; ) 



Procedure Genera-ePcp; 

begln ^c^z- := I -c PopSize do 

^ " f « ^cK^ - ,andc.(MIelDi.ension 

+ I ) ; 

end; > 

f 

•func --on MutaticnRate: word; 

ie ^ in Mutationals := Round(100-Sqri(I-?o ? Fitnsssi?cpCounirj)); 

end; > 

; 

Procedure Ai loca tePcpFi tness ; 

Var T emp_Fiiness : array [I . ■ PopMaxSize] of single; 

b£ - in for PopCountr :- 1 to PocSize do Teru P _Fitness f PcpCountr ] :- 0; 



for PccCcuntr := i to PopSize 



do 



for AiieiCountr := i" to GeneSize do 



Tem F itness PopCountr] T^Temp Fitness [ PopCountr] + 
k5S(Po ? List [PopCountr, AlleiCoun.r; - 
ModeIGene[AlIeICountr j ) ; 

for PopCountr := 1 to p ^pSize do Fitness [ PopCountr ] := 

if Temp_Fitness [PopCounir] - u men 

1 ro ^ r ^ nn1 -^i •= 0 9/Temo Fitness [ PopCountr ] 

else Temp_fitness [PopCountr J . u.s/i<=. 



ena; 
{ 



Function OneGeneFi tness : Boolean; 
TempOneGeneFitness : single; 



var 
begin 



f^lfe^rcu^tr 6 -^ l^GeneSize do 
f ° r X^SneGeneFitness := ^empOneGeneFxtness + 
A3S ( PopList { PopCountr, AllelCounu. ] 



Mod eiGene [ Al I un:r] } ; 

if PopFithess [ PopCountr j >= Terr.oCneGeneFit ness then OneGeneFi tn.ess 
: - False 

else OneGene Fitness := True; 

end; 

{ , 

Function AlielMutaticnRate : wore; 
begin 

if MutaticnSchema [AiieiCountr ] < = 0 then AiieiMutationRate : = 
AiieiMutationSeed 

else AiieiMutationRate := 

AiieiMutationSeed * (I + 
Trunc (Scrt ( Mut at icnSchema [ Ai ieiCount r ] } ) } ; 
end ; 

{ _ , 

Procedure PopMut a t ion-; 
var 

AiieiMutsticnDirection : integer; 

begin 

for PcpCcup.tr : = i to PopSize dc 

if Random (100) < MucaticnRate then 

for AiieiCountr := 1 to GeneSize dc 

if Random (100) < AiieiMutationRate then 
Begin 

PopList [PopCcuntr, AiieiCountr] : = 

Random ( Al lei Dimens-ion+i ) ; 

if CneGeneFitness then 
I no (Mutatiop.Scherua [AiieiCountr] ) 

else Dec (MutationSchema [AiieiCountr] ) ; 

end; 

end; 

; _____ _____ } 

Procedure CrossCver; 
var 

PopChi idren : word; 
CrcssOverLocat ion : word; 
GeneOne, GeneTwo : word; 

POpLack ; word; 

begin 

If Next PopCountr * 3 > 2 * PopSize then 

NextFopCountr (2 + PopSize) div 3; 

PopChildren := Next PopCountr + 1 ; 
PopCountr := 1; 

While PopCountr <= (NextPopCountr - 1 ) ' do 
begin 

CrossOverLocation : = Random (GeneSize - i) + 1; 

for AiieiCountr := 1 to CrossOverLocation do 
PopList [PopChildren, AiieiCountr] := 

PopList [PopCountr, AiieiCountr] ; 

Inc (PopCountr) ; 

For AiieiCountr := (CrossOverLocation + 1) to GeneSize 



do 



PopList [ PopChildren, AiieiCountr] := 

PopList [PopCountr, AiieiCountr] ; 




Inc (PopCountr) ; 

Inc (PopChilcren) ; 



for PopCountr :^ PopChi idren to PopSize do 
bee in 

GeneOne : = Random (Next PopCountr ) ^ I; 

GeneTwc : = GeneOne; 

While GeneTwc = GeneOne do 

Gene Two : = Ra ndom (Next? o pCou n zr) + 1 ; 

CrossOverLocat ion := Random ( GeneS i ze - Ij + 1; 

for AileiCour-ir : = I to CrossOverLocat icn do 
?oclist[?cpCountr, Ai lelCou.nc r ] := 

PopList [GeneOne, AiielCountr] ; 

For AlieiCcuntr := (CrossOverLocat ion 1) to GeneSi: 

Of f sprincList [ PcpCountr , AiieiCcuntr] : = 
PopList [ Gene Two , AiieiCountr] ; 

e n d ; 



end; 
{--- 



Procedure Selection; 
var 

AddGene : word; 

MaxGeneFitness : single; 

MinGeneFitness : single; 

RecrccucmcSpicies : array [ 1 PcpMaxSize] of word; 

TerTiO ?c o ' : ar ray [ 1 . . PopMaxSi ze , 1 . . GeneMaxSi ze ] of 



were; 
begin 

TotFitness := 0; 

Next PopCcunt r 
MaxGene Fitness 
MinGeneFitn 



0; 

PopFitness [ 1 ] ; 
PcpFitness [ 1 ] ; 

for PopCountr := lto Popsize do 
begin 

TotFitness :« TotFitness + PopFitness [ PopCountr ] ; 
if PopFitness [ PopCountr ] > MaxGeneFitness then 
Ma xGene Fitness : = PopFitness [ PopCountr] 

else if PopFitness [PopCountr ] < MinGeneFitness then 
MinGeneFitness : = PopFitness [ PopCountr] ; 

end; 

TotFitness := TotFitness / PcpSize; 

writeln(f generation, * , ' , MaxGeneFitness , ' , * , MinGeneFitness , 
ot Fitness ) ; 

if (MaxGeneFitness = 1) and (TotFitness >= 0.9) then 
begin 

Writ eOut ; 

Close (GenerationOutFile) ; 

Close (f); 

Halt; 



e rid ; 



enc ; 



for PopCountr := 1 to PopSize do 

if (MaxC-eneFitness - PcoFitness [ PopCountr ] ) < = 

(MaxGeneFitness - MinGene Fi t nss s } / Thres hoid2 Breed then 
begin 

Inc {Next PopCountr ) ; 
ReprocucincSpicies [ Next PopCountr ] := PopCountr; 
end; 

if NextPopCcuntr = 0 then 
begin 

writeln (' Population died out. Search not completed') 
Close (f J ; 

Close (Genera ticnOut File ) ; 
HALT; 

end; 



for PopCountr := 1 to Next PopCountr do 
for AiielCountr := 1 to GsneSize do 

Terr.pPop r PopCountr, AiielCountr] : = 

PcpList [ Reproducing Spicies [ PopCountr] , AllelCount 

for PopCountr := 1 to Next PopCountr do 
for AiielCountr := 1 to geneSize do 

PopList [PopCountr, AiielCountr] := 
TempPop [ PcpCour.tr, AiielCountr] ; 

if NextFcpCour.tr <> Trunc ( Next PopCountr /2 ) *2 then 
.begin 

Inc (Next: PopCountr) ; 

hddGene := Random ( Next PopCountr - 1) + I; 
for AiielCountr : = 1 to GeneSize do 

PopList [NextPcpCountr , AiielCountr] := 
PopList [AddGene, AiielCountr]; 

end; 



bee in 



MAIN ROUTINE 



Randomize; 

Assign ( f , * c : \algrthms\ga\ fitness . out * ) ; 
Rewrite (f ) ; 

PcpSize .:= 300; 
GeneSize := 10; 

InitiatelO; 

GenerateModelGene ; 

write (GenerationOutFile, 'Model Gene: '); 



fcr Alielcl^Pfr := i to Ger.eSiza do _ 

write (GenerationOutFiie, McdeiGene ( AlieiCoun t r } : 2) 
wr ite In (Genera tionOut File) ; 
writein (GenerationOutFiie) ; 

General ion := I; 

Genera te Pop; 



While Generation < MumOf Genera t ions do 
begin 

NuIlArrsys; 
AIlGcateFopFitness ; 

if Generation = 

Trunc ( Gene ration /Wr in e Frequency ) * Wri z e c rscuenc y 
then WriteOut; 

Selection; 
Crossover ; 
PcpMut anion ; 

writeln ( ' Generation : ' , Generation) ; 

if Generation = 

Trunc (Generation /Mutat ionSchemaBin) * 

Mu tcticnSchemaBin then 

for AlielCountr := 1 to GeneSize do 
. MutationSchema [AllelCountr ] 0; 

inc (generation) ; 

end; 



end . 



WriteOut; 

Close (Generat ionOut File ) ; 
Close ( f )*; 



APPENDIX 




ring . ini 

Apendix A and 3 Initialization file 



20 Number of APs 

QQ Lowest sensitivity for Rx in dBm 

2 Mininmal number of A? on ring 

4 Maximal Number of APs on the Ring 

2 upper bound of rings on A? 

2 1 - random position, 2 - grid posit 



ci : \sim\mac-sim\SGHz\Ring_dBm. T3L 
d : \sim\mac-sim\5GHz\Ring . cut 

i : \users\comnion\zvika\aviw\Ring_dEm. T3L 
i : \users\common\zvika\aviw\Ring . out 
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APPENDIX E 



Ring. out 

Transmitter to Subset Mapping Subset connectivity list 



Ring 


i : 


8 


3 


7 


9 


Ring 


2 : 


2 


3 


7 




Ring 


3 


1 


2 


8 


6 


Ring 


4 


10 


9 


6 


5 


Ring 


5 


12 


1 


10 


11 


Ring 


6 


14 


12 


5 


13 


Ring 


7 


15 


11 


13 


17 


?;ing 


3 


: 18 


17 


14 


19 





2 


3 


4 






1 


3 










2 


4 


5 




1 


3 


5 


6 




3 


4 


6 


7 




4 


5 


7 


3 




5 


5 


8 






6 


7 







Page 1 



APPENDIX F 



function cbwi=C3WL() 



# 



naces a channel reservoirs c 
p.ec by the example in Figure 



% T h i s fur.cz ion computes a 
": c 1 i cue £ . The cliques are 

5. 

7: unc:icn's input: 

'c-D Mapping vector between transmitters and Subsets 
Subset graph, defined as a connectivity matrix 
A Carrier-sense matrix that defines levels cf 
interference between the transmitters 
transmitter load vector - a matrix that defines the 
ccmmunicat icn intencity of each transmitter in terms 
cf frames per second 



12) 
^3> 

i 

I 
h 

'j. "unction 1 



s output 





Complete time/event r 


isp of the 




foil o i no e x amp 1 e : 




■5 


Time 


Tx-I Tx-2 


. Tx-21 


% 




i i 


i 


% 

•:■ 




i 6 


2 


:■ 




S 


S 










12) 


Each 


local reservoir 


behavior, 


s 


f o 1 1 owi nc e x amp 1 e : 




■V 












Time -## 




* 


Res-I 


Res-2 . . . Res -6 




A 


A A 




% 


3 


B 






C 


C 






D 


D 




% 









system, as presented in the 



as presented in the 



subnet s=zeros (6,21) / 

h Filling the subnet matrix 



with he relevant data 



n 


I, 


I) 


_ n . 
— / 








n 




2) 


_ o . 
— ^_ / 








n 


1, 


3) 


^ . 

— / 


n(2, 


3) 


= 1; 


n 




4 ) 


= 1 ; 


n(2, 


4 ) 


= 3; 


n 


1, 


5) 


= 3; 








n 


1, 


6) 


-3; 


n{2, 


6) 


= 5; 


n 


1, 


7 ) 


= 5; 








n 


1, 


8) 


-5; 








n 


I, 


9) 


= 5; 








n 


1, 


10) 


= 3; 


n{2, 


10) 


= 5; 


n 


i, 


11) 


-1; 


n{2, 


11) 


-3; 


n 


1, 


12) 


= 1; 


n<2, 


11) 


= 4 ; 


n 


1, 


13) 


= 4 ; 








n 


I, 


14) 


-4 ; 


n<2, 


14) 


«6; 


n 


1, 


15) 


= 6; 








n 


1, 


16) 


= 6; 








n 


I , 


17) 


-6; 








n 


1, 


18) 


= 4 ; 


n(2, 


18) 


= 6; 


n 


1, 


19) 


-1 ; 


n(2, 


19) 


= 4 ; 


n 


1, 


20) 


= 1; 


n(2, 


20) 


= 2; 


n 


1 , 


21) 


= 2; 









vNeigbor Clique connectivity graph 

ci ique= zeros ( 6 ) ; 

==1:6; 

ciiaue ( z , z ) =1 ; 
ciiaue (1, 2) =1; 
clique (2, 1) =1; 
clique (3, l)=i; 
ciiaue(4, 1)=1; 
clique(5, 3 ) = 1 ; 
ciique(6, 4)=±; 



clique (1, 3) =1; clique ( 1 , 4 ) =1 ; 

clique(3,5)=l; 
clique<4, 6)=1; 



. . / 

^Initial reservoir cf ag|^ciiques 
r.u;"_c:icue; = 6; ^For th^MJecific example 

-It is t:n.cc r t a n.t - o asSP less channels than transmitt 
' : - 1 . e . , otherwise nc reservoir nulling will occur 
r:um_channels = 10 ; 
reser voi r^cnes ( num_ci iques } ; 




^ As sign initial reservoirs 
reservoir {:} =num channels; 



^Channel request vector 
| channel_neec= zeros (1,21) ; % Tor the specific example case 

.-Channel request mienciiy expressed threshold 
channei_Neec_r hresholc= . 5; 

^Associate output files 

f 1- "c : \dat a \s incut \ reservoirs .cut"; 

f l="c: \data\simout \ channel_need . out " 

^Initial ti~e anf the total simulation duration in event units 

time=I; 

fin=5000; 

vMain loop. Perform execution for time=Q to time=fin 
} ■-, A z each point calculate the reservoirs of every clique 
"■and print it out 
for time=l : fin 

I Initiate the reservoir vectors" - it is assumed that each borrowing 
'^duration is equal to -one time step 
reservoir ( : ) =num_channei s ; 

Uirst, generate the channel request vector 
channel_".eed ( : ) =in ( random ( J , channel_Need_t hreshcid ) ; 

'^For ail nodes {transmitters in the network, chack whether they 
oneed a channel and then borrow one - iff available and update all 

r e s e r vo i r s accordingly 
for node= 1 : 2 1 

if c h a n n e l_ne ec (node) >0 

'VCheck to what cliques this node (transmitter) is associated 
t-and if possible borrow a channel 
if reservoir ( n ( 1 , node )) >0 

reservoir ( n ( 1 , node ) ^reservoir ( n ( 1 , node ) - 1 
else if n(2,node)>0 reservoir (n (1, node) ) >0 

reservoir (n ( 1 , node ) -reservoir ( n { 1 , node ) - 1 ; 
l end; 

end; % reservoir not emply 
pr intf ( fl , reservoir ) ; 
pr int f £ f 2 , channel_need) ; 
end; ?.channel_need 
end; for loop 
end; umain for loop 
close { f 1 ) ; 
close { f2 ) ; 



return; 
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